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EXECUTIVE  SUMMARY 

The  need  for  propellants  capable  of  delivering  Isp’s  of  600  seconds  or  better  has  led 
to  an  investigation  into  high  energy  density  molecules.  Such  species  may  derive  their 
energy  from  electronic  excitation  and  metastability,  unusual  bonding  situations,  or  from 
high  energy  ground  states.  The  Air  Force  is  currently  sponsoring  research  in  this  area 
through  the  AF  Astronautics  Laboratory,  the  AF  Office  of  Scientific  Research,  and  the 
AF  Aero  Propulsion  Laboratory.  A  Technical  Panel*,  funded  by  AFOSR,  administered  by 
Dr  Joseph  Early  of  Georgetown  University  and  chaired  by  Dr  William  Stwalley  from  the 
University  of  Iowa  provided  guidance  in  proposal  selection  and  insight  in  the  form  of  peer 
review  to  the  Air  Force  Steering  Group.  The  work  involves  both  theoretical  and 
experimental  investigations  and  some  of  the  nations  leading  research  teams.  In  order  to 
derive  synergistic  benefits  from  the  two  types  of  research,  the  first  annual  High  Energy 
Density  Matter  Contractors  Conference  was  held  in  May  1987. 

Technical  issues  facing  the  researchers  include  formation  mechanisms,  the  amount 
of  energy  stored  vs  the  molecular  weight  of  the  reaction  products,  duration  of  the  energy 
storage  (lifetime),  methods  of  stabilization  and  storage,  novel  bonding  structures,  and 
decay  mechanisms.  The  work  can  be  loosely  grouped  under  the  following  categories: 
screening;  bonding  and  reaction  studies;  theory  verification  and  specie  characterization; 
and  feasibility/stability  demonstration.  Theorists  are  performing  quantum  mechanical 
calculations  to  determine  energy  states  and  predict  possible  decomposition  pathways. 
Experiments  are  underway  to  verify  theory,  characterize  species,  and  investigate 
stability.  Several  of  the  important  conclusions  presented  at  this  meeting  are  reported 
here. 

An  energetic  stable  state  of  H4  had  been  predicted  by  Nicolaides  in  his  maximum 
ionicity  excited  states  (MIES)  theory.  Initial  studies  showed  that  an  energy  minimum 
existed  in  the  ground  state  and  first  excited  state  of  the  proposed  molecule.  After 
extensive  investigation,  the  ground  state  of  tetrahydrogen  (H4)  was  found  to  be  an 
unbound  saddle  point.  This  was  confirmed  independently  by  the  groups  from  UC  Berkeley 
(Lester),  Unitdd  Technologies  Research  Center  (Michels),  and  the  US  Army  Ballistic 
Research  Laboratory  (Lengsfeld).  Two  other  MIES  candidates,  HLi3  and  the  excited  state 
of  H4,  still  show  promise  of  stability  and  high  energy  content. 

Several  groups  reported  on  powerful  new  computational  methods  developed  to 
investigate  energy  configurations.  Their  almost  simultaneous  independent  findings  on  H4 
lend  great  credence  to  the  concept  of  "theoretical  experiments"  being  done  on 
supercomputers.  The  new  computational  tools  allow  investigators  to  quickly  and 


accurately  pred'ct  the  potential  stability  of  a  candidate  system  before  it  is  prepared  in 
the  laboratory.  Once  a  stable  state  is  identified,  formation  and  decomposition  pathways 
can  also  be  calculated.  This  greatly  reduces  the  time  and  risk  involved  with  laboratory 
experiments  on  the  species. 

Michels  reported  calculations  on  a  new  energetic  oxidizer  predicted  to  be  stable  in 
the  ground  state.  This  compound,  asymmetric  dinitrogen  dioxide  (a-N202),  is  an  analogue 

*  of  fluorine  azide  (FN3)  and  is  predicted  to  be  a  cryogenic  solid. 

Benard  reported  that  FN3  can  be  viewed  as  a  loosely  bound  complex  of  N2  and  the 
metastable  excited  singlet  NF.  He  found  that  by  using  FN3  as  a  starting  material,  one 

#  can  generate  high  concentrations  of  the  excited  singlet  NF,  a  short  wavelength  laser 
species. 

Yarkony  presented  a  recently  developed  algorithm  for  treating  spin  interactions  in 
the  computation  of  radiative  lifetimes.  This  method  enables  previously  untenable 
calculations  to  be  performed  which  help  predict  the  lifetimes  of  candidate  species. 

Helm  reported  on  a  newly  developed  technique  for  experimentally  observing  high 
Rydberg  states.  This  rather  elegant  method  is  used  to  precisely  characterize  the  energy 
surfaces  of  candidate  molecules. 

These  items  represent  some  of  the  positive  results  gained  in  the  first  year  of  work 
by  the  initial  HEDM  investigators.  There  are  now  over  twice  as  many  efforts  underway, 
expanding  our  view  into  new  theoretical  and  experimental  approaches.  While  our  main 
objective  is  to  identify  and  develop  new  high  energy  density  molecules,  it  is  important  to 
note  that  we  are  also  advancing  the  state  of  the  art  in  chemistry  and  chemical  physics. 

♦The  Steering  Group  has  established  a  Technical  Panel  to  consult  on  the  scientific 
merit  of  proposals  and  to  provide  peer  group  review  of  technical  progress  as  research 
investigations  proceed.  Dr  William  Stwalley,  University  of  Iowa,  a  noted  specialist  in  both 
chemistry  and  physics  chairs  the  Panel.  Other  Panel  members  are: 

f 

Dr  Charles  F.  Bender  Advanced  Computational  Methods  Center,  Athens,  GA 

Dr  William  Hopper  Princeton  University 

ft 

Dr  M.  Frederick  Hawthorne  UCLA 

Dr  Ron  R.  Herm  Aerospace  Corporation,  Los  Angeles,  CA 

Dr  Marilyn  Jacox  National  Bureau  of  Standards 

Dr  Lewis  H.  Nosanow  National  Science  Foundation 

Harvard  University 


Dr  Isaac  Silvera 
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~  MATRICES  ARE  HELD  TOGETHER  BY  RELATIVELY  WEAK 
DISPERSION  FORCES 


Specific  impulse 


UNH 


«U  P  u  w 


iiimywtirn 


BTunwnieiHWVfw r*  irr  rww^onr* 


A 


% 


(0 

Q. 

E 

o 


>< 

cr 

^  h- 
<0  < 
X  5 

«  LU 

H 

1  m 

«  u. 
o>  o 
<o  w 

^  Z 

+  o 

I  o 
I  < 
^  oc 

CO  li 

X  . 


rs.  o 

co  5 


rr  « 
K  O 

<  m* 
s  2 
S  ° 


3  2 

3  s 

5  < 

i  = 

o  g 

iu  __J 


©  co 
«o  a> 

O  N 


O  N 

•  • 

O  40 


in  i^» 
o  n 

©  in 


cm  oo 

o  co 
•  • 

o  co 


I  I 


X  X 

<3 


Gaseous  Energy  Content 


14 


UVUWffUWVUVUV 


WUVVVJI UVUVUIUVIIVUVI 


CO 

0) 

3  O 

(0  w 

CO 


O  s 


< 

o 


c 

0) 


h 

CO 

Ui 


O  jjj 

C/>  m 


>.  o 

CD  •"" 

* 


p» 

(0 

UJ 


x 

b- 

< 


CO 


o 

h- 

< 

LL 

o 


CO 

UJ 


< 

X 

t- 


LU 

UJ 


UJ 

x 

< 

i- 

< 

x 

£ 


c^* 

C5 


O  UJ 


CO 


C“ 

UJ 

o 


o 

LL 

o 


UJ 

D 


CO  ... 

m  UJ 

O 


o 

X 

UJ 


UJ 

X 

LLi 

a 


ill 

a 


UJ 

X 


UJ  UJ 
CL 


x  r- 


O 

P 

< 


UJ 

Q 


fr“ 

O 

< 

UJ 


UJ 

N 

CO 


p  UJ 

t  ° 

<  :d 

£  t 

i  2 

UI  <[ 


H-  H* 


UJ 

X 

< 

K 


X 

X 

X 

UJ 

UJ 

UJ 

X 

X 

X 

h- 

H- 

h 

22 

CO 

co 

»- 

H 

H* 

< 

< 

< 

X 

X 

X 

£ 

£ 

£ 

i 

1 

1 

15  '16 


>- 

< 


x 

x 

X 

h- 

< 


CO 


o 

LL 


ill  o 


h- 

< 

X 


LU 

CL 

CO 


Z  UJ 
2  £ 


Z  <  CO 


N  UJ 

X  o 


UJ 


o 

UJ 


u_  & 

It  UI 


UJ 

X 

< 


o 

xO 
<  c 
x* 

i_  UJ 

<co 


UJ 


£  Si= 


UJ 

X 

f- 


Q 

UJ 

NJ 

5 

ca 

< 

h- 

CO 


>CO 
Z  UI 
<  CO 


UJ 

ffi  _ 

CO  UJ 
UJ  > 

—  UJ 

Oil 

UJ 

X  x 
CO  UJ 


o  s 

—  Q 


oS 

I-  o 

<  CC 

UJ  UJ 


< 

UJ 

£ 


o 


CO 

u. 

o 

CO 

o 


CO 


X 

UJ 


h 


o 

< 

<« 

oS 


ui  o 
?x 


UJ  o 
X  o 

<  i 

Fo 
>  m 


<o 

OS 


4 


§ 


3 


Towards  Understanding  the  Stability  of  the  Cluster 

A.  Netropoulos  and  C.A.  Nicolaides 

Theoretical  and  Physical  Chemistry  Institute, 

National  Hellenic  Research  Foundation,  Athens,  Greece 


As  a  first  attempt  to  study  the  stability  of  the  H,*(C3V)  cluster 
we  have  used  the  MRD-CI  methods  and  a  medium  size  basis  set  to 
calculate  various  sections  of  the  potential  energy  surfaces  of  its 
ground.and  first  excited  states.  We  show  that  these  correlate  to 
the  (XXE  TI  +)  and  (X1!  Vi  +)  states  of  the  two  H2 
constituents  respectively?  Finally,  we  report  on  the  calculation 
of  the  diabatic  matrix  elements  of  the  vibronic  interaction  in  the 
region  of  the  avoided  crossing  which  is  crucial  to  the  stability 
of  H4*. 


ON  THE  STABILITY  OF  THE  H4  EXCIMER 

A.  Metropoulos  and  C.A.  Nicolaides 

Theoretical  and  Physical  Chemistry  Institute 

National  Hellenic  Research  Foundation 

48  Vas.  Constantinou  Ave.,  Athens  11635,  Greece 


HISTORICAL  PERSPECTIVE 

In  our  institute  there  has  been  a  long  standing  emphasis  on 
the  study  of  atomic  and  molecular  excited  states.  Within  the 
course  of  these  studies,  it  was  decided  in  1982  to  compute  the 
low  lying  excited  states  of  the  He-H,  with  the  objective  of 

«*•  A 

doing  dynamics  on  these  surfaces.  The  idea  was  to  do  more 
accurate  calculations  than  those  reported  in  the  literature  /I/ 
and  to  map  a  wider  portion  of  the  hypersurface  including  the  H-H 

Q 

stretch  and  a  45  section.  What  was  expected  was  tne 
construction  of  accurate  ground  and  excited  states  van  der  Waals 
surfaces  for  energy  transfer  calculations.  In  the  course  of  the 

investigation  it  became  evident  that  for  a  long  enough  H-H 

o 

distance  (4.0bohr)  and  a  45  angle  there  was  a  non-van  der  Waais 
minimum  of  about  1.52  eV  with  respect  to  He  +  H^.  The 
corresponding  He-H^  distance  was  about  1.5  bo'nr,  or,  else,  the 

closest  He-H  distance  was  1.41  bohr.  In  this  stable  geometry 

>--**■ 

there  is  an  avoided  crossing  between  the  He  +  H^(3  4-w)  excited  and 
the  He+H^CX  1^)  ground  states  as  the  H-H  distance  varies  / 2 /  (not 
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I 

I 

•f 


present  in  colAinear  or  geometries).  Figure  1  shows  the 
potentials  of  the  ground  and  the  first  excited  states  near  the 
minimum.  The  existence  of  such  an  avoided  crossing  has  explained 
for  the  first  time  the  observed  fluorescence  quenching  of  HD  / 3 / 


very  satisfactorily.  The  discovery  of  a  geometry  for  HeH^  with  a 


chemical  minimum  gave  rise  to  the  so  called  maximum  ionicity 
excited  state  (MIES)  theory  as  a  means  of  explaining  the 
existence  of  such  a  minimum  in  terms  of  general  chemical 
features.  The  impetus  for  such  a  theory  was  the  fact  that  the 
H,(B  2-h.)  state  is  ionic  at  intermediate  H-H  distances  with  the 


maximum  ionicity  occuring  at  1.0  bohr/i/.  Thus ,  as  He  approaches 


H-  the  H-H  bond  stretches  close  to  the  maximum  ionicity  iinil 


which  favors  energy  reduction  through  small  overlap  between  na~ 
and  He  and  a  large  Coulomb  interaction  between  H^*"  and  He. 


Support  for  this  model  comes  from  the  fact  that  the  equilibrium 


distance  of  the  HeH  molecular  ion  is  1.11  bohr/15/  which  is  very 
close  to  the  closest  He-H  distance  at  the  equilibrium  geometry. 


Extending  this  model  of  bonding  to  larger  systems 


interacting  singlet,  closed  shell  moieties,  one  can  make  very 

t  h  e  t.  r 

useful  predictions  of  approximate  "stable"  geometries  of, excited 


surfaces.  Taking  H^  as  an  example,  this  scheme  can  be  described 


as  follows.  An  excited  state  of  clusters,  say 


ABCD )  , 


correlating  with  the  fragments  ( AS )+( CD )  ,  is  thought  of  as 

S*-  J~- 


occuring  by  an  interaction  of  the  tyre  (A! 


The 


geometry  of  the  bound  (ABCD)  complex  is  obtained  from  the 


geometry  of  (ABC)  ,  and  by  placing  D  a  distance  from  C  at  wh 
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to  the  wave 


M 


ionic  V3  structures  have  a  maximum  contribution 


function  of  the  (CD)  fragment.  For  H^,  the  geometry  of  the 


positive  ion,  H*,  is  an  equilateral  triangle  with  side  of  1.65 
bohr  / 5 / .  The  fourth  H  atom  is  then  placed  above  the  centre  of 


the  triangle  at  a  distance  of  about  4.0  Bohr,  which  is  the 
maximum  ioncity  distance  of  the  B  Z-u.  state  of  H,  / 4/ -  Figure  2 


shows  the  relevant  geometry  and  the  coordinate  system  used.  The 


MIES  prediction  on  was  supported  by  the  results  of  MRD-CI 


calculations,  which  gave  a  minimum  at  the  triangle  side  of  1.70 
bohr  and  a  distance  of  3*8  bohr  from  the  centre  of  the  triangle 
for  the  fourth  HI  /6/.  The  calculations  were  done  in  Cg  symmetry 
and  figure  3  shows  the  potentials  of  the  two  lowest  A'  states  f:r 
various  lengths  of  the  side  of  the  equilateral  triangular  base. 
Notice  that  while  work  on  the  excimer  has  been  done  by  other 
workers  for  various  geometries  /7/,  the  application  of  the  MIES 
scheme  has  lead  directly  to  a  pyramidal  geometry  as  the  minimum 


energy  geometry  for  the  excimer.  All  previous  considered 


geometries  lead  only  to  saddle  points  of  the  K  hypersurface  / 8 / . 


FURTHER  INVESTIGATION  OF  THE  K.  EXCIMER 

4 


After  the  MIES  theory  had  been  succesfully  tested  / 6/ ,  it 

was  decided  to  investigate  more  accurately  the  H  .  excimer  state 

4 

because  we  needed  its  character  and  its  accurate  wavefur.ctions 
and  surfaces.  As  previously,  the  MRD-CI  method  /9/  was  employed 
and  the  calculations  were  done  in  C^.  symmetry,  due  to  program 


20 


VI 


•I 


*>.  *1.1^  *  m.  r  . 


limitations,  but  now  a 


•ger  tasis  set  wa. 


/TO/. 


attention  was  focused  on  the  following:  a)  Obtain  more  accurate 


surfaces  ana  wav 


e  f  u  n  c  t i o  n  s  while  maintaining  the  C 


3  V  3ec~e -ry . 


b)  Use  these  wavef  unctions  tc  obtain  the  r.cr.-aiiabati  o  matrix 


elements  of  the  vibrcr.ic  interaction  in  the  vicinity  of  th; 


avoided  crossing,  assuming  all  motions  frozen  except  the  cr.e 

along  R  (see  fig.  2).  c)  Verify  that  the  first  two  *A'  states  of 

4-  j-  +  1  y  r  i  r  f  1  t 

H  correlate  to  the  (X  g  X  ^—g)  and  (X  i-  g  3  2-  u)  states 


of  two  hydrogen  molecules 


At  this  point,  the  approximate  rature  of  the  MIES  bases  pr- - 


diction  of  the  geometry  of  the  excited  state  minimum  became  evi¬ 


dent.  As  it  was  pointed  out  to  us  /I  4/  and  a; 


ran  be  deduced 


from  the  calculations  in  ref.  8,  if  the  true  C?,.  symmetry  i. 


considered,  there  is  no  avoided  crossing  between  the  first 


totally  symmetric  states  but  a  true  crossing  between  an  A 


an  E  state  (E  =  A'  +  A").  The  geometry  of  tr.is  crossing  was 


later  calculated  in  C ^  symmetry  /'2/  and  it  was  found  to  be 

1.8  bohr  and  R  =  4.225  bchr  at  an  energy  cf  -2. 1617  a.u.  Fi 

i  d. 

4  shows  the  extracolateu  curves  of  the  A,  ana  a  states  tr.e 


1 


intersection  of  which  has  been  taken  as  the  crossing  ccir.t. 


viousiy  then,  the  stable  geometry  must  be  of  a  tr\ 


From  the  optimization  calculations  with  respect  tc  motions  of  H, 


in  C3y  geometry,  represented  i 


n  figure  5,  it  was  deduct 


Zg  symmetry  has  to  be  arrived  at 
angular  base  rather  than  by  disc 


searching  for  the  cptim..m  g*  of  *_• 


£**»Tr*' 


Here,  we  consider  only  the  deformation  of  the  equilateral 
triangular  base  into  an  isoskeles  (but  not  into  a  scalene)  tri¬ 
angle  /12/.  After  a  sequence  of  optimizations,  we  arrived  at  a 
geometry  having  a  true  avoided  crossing  between  the  X  AA'  and  the 
A  A*  states  which  is  about  154  meV  lower  than  the  (  A^,  E) 
crossing  point  (-2.0674  a.u.  for  the  A  *A'  and  -2.0707  a.u.  for 
the  X  *A'  (state)  and  which  is  a  minimum  of  the  A  ^A*  state.  At 
this  geometry,  the  equal  sides  of  the  isoskeles  triangle  are  1 .78 
bohr  each,  the  third  side  is  1.62  bohr  and  is  situated  on  the 
z  axis,  3.90  bohr  above  the  plane  of  the  triangle.  The  origin  of 
the  coordinate  system  is  taken  at  a  point  2/3  from  the  vertex  of 
the  perpedicular  bisector  of  the  isoskeles  triangular  base.  Figu¬ 
re  6  shows  the  variation  of  the  potential  with  R  for  the  opti¬ 
mized  triangular  base. 


Regarding  the  correlations 
shows  that  they  are  correct  as 
symmetry  in  this  table  is  Cjy  , 


previously  mentioned,  table  1 
stated.  Although  the  final 
the  correlations  are  valid  for 


the  optimized  C^.  symmetry  as  well.  This  is  so  because  in  this 


case,  the  A^  irrep  of  correlates  uniquely  to  the  A’  irrep  of 

C  =  &~v)  /I  0/ . 


Finally,  the  non-adiabatic  matrix  elements  for  the  radial 
coupling  were  computed  in  the  neighborhood  of  the  avoided  cross¬ 
ing  for  the  optimized  H  *  (C^.  )  geometry.  The  finite  differences 
method  of  Lorqueo  and  coworkers  as  it  has  been  modified  for  the 


except  for  the  r  internal  coordinate .  Table  2  gives  tr.e  result 
of  this  calculation,  "coice  that  while  the  strer.gtr.  of  the  inter¬ 
action  is  snail,  it  persists  for  a  rather  large  interval.  The 
usefulness  of  this  result  is  that  one  can  sake  ar.  order  of  magni¬ 
tude  estimate  o:  the  transition  rate  to  the  ground  dissociative 
state  by  using  methods  already  developed  for  diatoms  / 13/. 


LIFETIME  OF  THE  H.  EXCIMER 

4 


A  commutation  of  the  lifetime  of  H  is  verv  desirable.  There 

V 

appear  to  be  three  mechanisms  by  which  K  may  dissociate:  a)  By 
a  radiative  transition  to  the  ground  dissociative  state.  The 
probability  of  this  transition  should  be  rather  small,  due  to  the 
proximity  of  the  levels  at  the  avoided  crossing  ar.d  due  to  the 
smallness  of  the  Frar.k-Condon  factors  away  from  it.  b)  Ey 
exciting  the  system  to  the  level  of  the  (  A^,  ^  E )  crossing  point. 
Thus,  the  system  car.  be  "stable’1  only  at  low  vibrational  states 
and  the  highest  such  state  must  be  found.  c)  By  a  diabatio 
crossing  to  the  ground  state.  The  probability  of  this  transition 
may  very  well  be  the  factor  cor.troling  the  stability  of  H ^  . 
Calculations  for  the  elucidation  of  the  first  two  mechanisms  are 
now  in  progress.  Crier  of  magnitude  calculations  regarding  the 
third  mechanism  (for  the  R  internal  coordinate  only)  are 
contemdated  for  the  near  future. 


rfy  *  iWWIrtr®? 
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Iac-e  1 


ce tween  * j 


the  suns 


lal 


energies 


r 


.  Comparison  between  1;  the  sums  of  the  potent 
of  two  hydrogen  molecules  at  infinite  secaraticr  and  a  T  shared  C-  , 

A  V 

complex  at  large  internolecular  separations  (1st  three 
columns)  and  2)  the  potential  energies  of  computed  with 

reference  configurations  initially  corresponding  to  K  -K~  ( C_  )  and 

Cs)  (last  three  columns).  The  MR  DC  I  potential  energies  of 

X 

the  molecules  are  also  given.  For  the  Cs  symmetries,  should 

be  reolaced  with  A':  &  is  the  angle  between  K  and  the  z  axis. 

^  is  the  bond  length  of  an  H  diatom. 


^  (a.u. ) 

H  ( r  = 

1.8,  R=4 

.2  a.u.) 

Geometry 

1  .4 

1.5  2.43 

Cs* 

Cs-f 

x'ZriZs 

-2.3400 

-2.3375  -2.1942 

-2.1418 

-2.1139 

-2-06:9 

X  1Ai 

-2.3393 

-2.3369  -2.1947 

-2.1421 

-2.1137 

-2.0620 

Ir+  1T  + 

X  Z.  g+3  d-u 

-1 .8735 

-1.8867  -1.8513 

-1 .9875 

-2.CC84 

—  •  VJ  C  •  *4 

A  ^A^ 

-1 .8730 

-1.836 3  -1.8506 

-1 .9852 

-2.CC80 

'1  o  '  -  ; 
-4.00 ■ ^ 

K  (X1?g+) 

^  1T'* 

H ^  (B  4u) 

-1 .16999 

-0.70350 

-1 . 16874  -1 .0971  1 

-0.71799  -0.75424 

*  ^  =  45  deg 

(Cs) 

+  ^  =  0  des:  ( C  ) 
jv 
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Table  2.  Non-adiabatic  matrix  elements  (in  a.u.)  of.  the  vibronic 
coupling  along  the  R  coodinate  (R  in  bohr). 
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Figure  1  .  The  potentials  of  the  ground  and  the  first 

as  functions  of  the  H-K  distance.  -  R  =  2.0  bonr; 

bohr . 


Figure  6.  The  variation  of  the  potential  with  R  for  the  X  anc 

the  A  "'a.  states  of  H,  (C  );  -  eigenvalues; 

-  extrapolated  energies.  The  horizontal  lines  show 

the  corresponding  levels  of  the  (1A1,"'s)  crossing 
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ABSTRACT 

These  studies  consist  in  attempting  to  detect  the  existence  of  a  metastable  1  A’ 
state  of  the  H4  molecule  and  in  determining  its  lifetime.  The  H4  preparation 
methodology  consists  in  crossing  an  intense  beam  of  metastable  H3  molecules  in  the 
2pz  2A"2  Rydberg  state  with  a  beam  of  ground  state  X  ^  HI  molecules.  The 
calculated  geometry  of  H4  (1  A')  is  a  pyramid  whose  base  has  a  nuclear  configuration 
close  to  that  of  the  2pz  2a*2  H3  molecules.  Furthermore,  the  height  of  that  pyramid 
is  approximately  the  same  as  the  effective  length  of  the  2 pz  Rydberg  orbital.  The 
reaction  H3(2pz2A"2)  +  HI  (Xi£*)  ■*  H4  PA’)  +  I  (2P3/2)  is  endothermic  by  about  1  eV 
and  the  H3  (2pz  2a*2)  beam  has  an  average  translational  energy  of  about  8  eV. 

The  crossed  beam  apparatus  has  been  assembled  and  put  into  operation.  The 
primary  H3  (2pz  2A”2)  beam  is  generated  in  a  high  temperature  arc  discharge  source 
through  which  H2  flows.  The  discharge  plasma  emerges  from  a  cathode  nozzle 
orifice,  and  the  beam  passes  through  a  skimmer,  then  through  a  differentially 
pumped  chamber  and  into  the  main  scattering  chamber.  The  secondary  beam  of  Ar 
(in  control  experiments)  or  HI  (or  Dl)  is  generated  by  expansion  of  the 
corresponding  gas  through  a  capillary  array.  The  scattered  species  are  detected  by  a 
rotatable  quadrupole  mass  spectrometer.  A  report  will  be  given  on  the  progress 
achieved  in  these  experiments  so  far. 
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.  UP  -  amusion  pumps,  PU  -  primary  beam  source 
the  plane  of  the  drawing),  Nz  -  nozzle,  SK  - 
skimmer,  VS  -  velocity  selector,  SB  =  secondary  beam  source  (with  the 
beam  perpendicular  to  the  plane  of  the  drawing),  CC  =  correlation 
chopper  blade,  MS  =  mass  spectrometer  detector,  TMBb  =  501/s 
turbomolecular  pump  for  buffer  chamber,  I  =  ionizer  and  ion-focus 
lenses,  CEM  =  channeltron  electron  multiplier,  CP  =  350 1/s  liquid  He 
cryopump,TMP  =  360 1/s  turbomolecular  pump,  IG  =  ionization  gauge. 
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Side  view  of  hybrid  apparatus  with  the  metastable  Ha  arc  discharge 
'primary  beam  source  connected  to  the  crossed  molecular  beam 
Apparatus.  vOl.  vacuum  chamber:  S.  water-cooled  copper  skimmer;  GV, 

fmeumatic  gate  valve;  F,  flexible  flange;  B5, 50-inch  bell  jar.  Dashed 
ines  indicate  box  within  the  bell  jar  which  contains  the  magnetic 
velocity  selector. 
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ABSTRACT 


The  2pz  2a*2  Rydberg  excited  state  of  H3  is  known  experimentally  to  be 
metastable,  but  its  predissociative  lifetime  has  neither  been  measured  nor 
calculated,  although  a  lower  bound  of  40  psec  has  been  experimentally  estimated 
for  it.  Furthermore,  the  potential  energy  surface  of  the  1A'  state  of  H4  has  been 
calculated  to  have  a  well  whose  depth  is  1.86  eV.  The  objective  of  our  theoretical 
studies  is  to  perform  calculations  of  the  predissociative  lifetimes  of  both  of  these 
molecules. 

The  initial  calculations  are  being  performed  for  the  H3  system  We  are  using 
initially  a  double  many-body  expansion  (DMBE)  of  the  two  lowest  potential  energy 
surfaces  of  H3.  The  computation  involves  three  steps:  a)  Scattering  calculations  on 
the  lowest  state.  We  have  performed  these  accurately  for  energies  up  to  16  eV 
using  hyperspherica!  coordinates  and  are  investigating  approximate  methods  for 
extending  them  to  the  higher  energies  required,  b)  Bound  state  calculations  on  the 
upper  surface.  The  structure  of  that  surface  is  being  analyzed  and  a  selection  is 
being  made  of  the  basis  functions  to  be  used,  c)  Interaction  matrix  calculations. 
These  have  not  yet  been  initiated.  The  results  obtained  so  far  will  be  presented  and 
discussed. 
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The  four-center  H4  system  has  been  the  subject  of  theoretical  studies 
since  the  late  sixties,  primarily  due  to  interest  in  the  relationship  of  the 
ground  state  hypersurface  to  the  H2  +  D2  *  2HD  exchange  reaction . 1 » 2 > 8 
Investigations  into  the  nature  of  the  excited  state(s)  of  such  system  have 
been  limited  in  number  even  though  SCF-MO-CI  calculations  in  1976  by  Goddard 
and  Csizmadia  indicated  the  existence  of  a  bound  "excimer  state. Recently, 
however,  Nicolaides  et  al. ,  using  an  extension  of  their  modeling  work  on  rare 
gas  hydrides8  have  predicted  the  existence  of  a  bound  singlet  excited  state  of 
H4  and  have  extended  this  analogy  to  the  (H2)2  and  (83)5  clusters.6*7 

The  theoretical  studies  of  Refs.  6  and  7  predict  that  the  bound  state  of 
H4  arises  from  the  reaction  of  molecular  hydrogen  in  its  ground  state  with  H2 
that  has  been  excited  to  the  state. 

H2(B,Eu+)  +  H2(X1Eg+)  +  M  -*•  H4*  +  M  (1) 

Further,  the  stable  excited  state  is  predicted  to  be  produced  only  when 
the  B1£u+  state  precursor  has  sufficient  "ionic"  character  for  reaction  with 
H2(X).  This  is  predicted  to  occur  for  v  >  3  in  the  upper  state.6*7 

Quenching  studies  of  HD(B,v)  performed  by  Atkins,  Moore,  and 
coworkers8-10  show  an  unusually  large  quenching  cross  section  (79A2)  for  H2. 
This  has  led  to  further  speculation  that  formation  of  a  complex  is  involved. 
The  H2(B)  state  in  v'  -  3  has  sufficient  energy  to  dissociate  itself  as  well 
as  an  H2  collision  partner: 

H2(B)  +  H2  +  4H  AE  -  -2.70  ev 

Hence,  quenching  of  H2(B)  by  H2  may  be  highly  reactive  with  H-atoms  as  the 
products . 

The  goal  of  our  experimental  program  is  to  produce  H2(B,v)  in  a  super¬ 
sonic  jet  by  several  techniques,  provide  a  sufficient  number  of  collisions 
with  H2  to  produce  H4* ,  stabilize  the  H4* ,  and  probe  for  its  existence  and 
molecular  properties.  Four  different  methods  have  been  chosen  to  produce 
H2(B,v)  in  the  jet:  1)  energy  transfer  from  Ar*  resonance  states  (1Q48A, 

1066A  lines)  in  Ar/H2  mixtures  via  discharge  production  of  Ar*  in  the  nozzle 
expansion  region;  2)  energy  transfer  from  Ar*  resonance  states  in  Ar/H2  mix¬ 
tures  via  resonant  excitation  of  Ar*  with  a  VUV  discharge  lamp;  3)  direct 
three  photon  laser  excitation  of  H2(b,v)  at  the  nozzle  throat;  and  4)  direct 
single  photon  laser  excitation  of  H2(B)  via  generation  of  tunable  VUV  radia¬ 
tion  from  third  harmonic  generation  in  rare  gases. 
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COLLISION  FREQUENCY  (s'l) 


The  primary  reason  for  using  a  supersonic  molecular  beam  to  produce  H4* 
is  to  provide  a  locally  high  density  but  optically  thin  source  of  H2  for 
direct  VUV  excitation  of  the  H2(B)  state  as  well  as  collisional  production  and 
stabilization  of  H4*  .  Three  primary  criteria  must  be  met  by  the  source: 

1.  There  must  be  sufficient  collisions  in  the  expansion  beyond  the  nozzle 
throat  for  collisional  quenching  of  H2(B)  to  occur  at  a  rate  at  least 
comparable  to  that  for  spontaneous  radiation. 

2.  When  energy  transfer  from  Ar*  to  H2  is  employed  there  must  be  sufficient 
collisions  beyond  the  nozzle  throat  for  Ar*  to  be  quenched  by  H2. 

3.  Formation  of  stabilized  H4*  is  a  three-body  process.  Hence,  sufficient 
collisions  must  occur  in  the  excitation  region  between  the  colliding 
H2(B)  -  H2  complex  and  a  third  body  (h2  or  Ar) . 

An  analysis  of  the  nozzle  flow  characteristics  (Figure  1)  show  that  H2(B) 
must  be  excited  within  0.9  mm  of  the  nozzle  throat  in  order  for  conditions 
1  to  3  to  be  satisfied. 


EXPANSION  01  STANCE  (roil 

A-4608 

Figure  1.  H2(B)  Collision  Frequency 

as  a  Function  of  Distance 
from  the  Nozzle  Throat 


The  calculations  of  Nicolaides 
and  coworkers5"*7  show  that  H4*  is 
bound  by  3.1  ev  with  respect  to  dis¬ 
sociation  to  form  H2(B)  and  H2 .  How¬ 
ever,  they  also  have  shown  that  the 
complex  is  bound  by  1 .86  ev  with 
respect  to  dissociation  along  the  H“ 

+  Hj"1"  coordinate.  Hence,  photolysis 
of  H4*  may  produce  H3+  +  H“ ,  H3  +  H, 
or  H2(B)  +  H2.  A  careful  selection 
of  photolysis  wavelengths  may  result 
in  selective  excitation  into  any  of 
these  channels.  In  our  experimental 
effort  we  will  be  using  a  range  of 
photolysis  wavelengths  in  order  to 
excite  the  channel  leading  to 
H2(B)  +  H2 .  We  are  also  considering 
detection  schemes  to  observe  ionic 
dissociation  pathways  and  H2(6) 
quenching  collisions  which  produce 
atomic  hydrogen.  Observation  of  H4* 
may  also  be  accomplished  via  detec¬ 
tion  of  passive  emission  from  the 
radiating  excimer,  if  such  radiation 
does  indeed  exist. 

Photolysis  of  H4*  win  initially 
be  attempted  using  a  broadband  flash- 
lamp  source  (X  =  190  to  3000  nm).  A 
solar  blind  PMT  will  be  used  to 
detect  H2(B)  emission 


resulting  from  the  photolysis.  Use  of  this  broadband  source  is  desirable  in 
order  to  minimize  the  time  required  for  a  long  search  for  appropriate 
photolysis  wavelengths.  As  H2(B)  emission  is  observed  the  flashlamps  will  be 
filtered  to  narrow  the  effective  photolysis  bandwidths  until  a  spectral  range 
is  defined  over  which  more  selective  laser-based  photolysis  studies  can  be 
conducted. 
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CHARACTERIZATION  OF  TETRAHYDROGEN  VIA 
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TETRAHYDROGEN:  NICOLAIDES  et  al . 


STABLE  EXCIMER  OF  H2(B)  +  H2 

FORMED  WHEN  H2(B)  IN  STATE  OF  "MAXIMUM  IONICITY" 
CORRESPONDING  TO  v‘  =  3  OF  B-STATE 

STABILITY  OF  H4 

H4  *  H2(B)  +  H2  AE  =  -2.48  eV 

H4  ♦  H3+  +  H"  AE  =  -1.86  eV 

VIBRATIONAL  FREQUENCIES  PROBABLY  SIMILAR  TO  H3+ 
SYMMETRIC  STRETCH:  v  =  3185  cm-1 
BEND:  v  =  2516  cm"1 


/  I  \ 


GEOMETRY 


EXPERIMENTAL  APPROACH 


THREE  PHOTON  EXCITATION  OF  H2(B) 


QUENCHING  IN  H2  AND  H2/Ar 


#  H2(B)  quenching  by  h2 

MOORE  et  al .  cth2  =  79A  2 


c  =  fiii) 

L  ^  Ilm  J 


1/2 


1/2 

^  s  1.1  x  10  ^  cm3  molecule  ^  s  1 

QUENCHING  CONDITION.  kq[H2]  =  krdd 
krad  =  1.5  x  109  s"1 

[H2]  =  1.4  x  1018  molecule  cm"8  AT  300  K 
•  A r*  QUENCHING  BY  H2 

BOXALL  et  al.  a  =  48A2 


-If)  t  3  -1-1 

kg  -  1.9  x  10  (-^qq)  cm  molecule  s 

QUENCH  A r*  1066A  LINE:  krad  =  1.2  x  108  s'1 

[H2J  =  6.3  x  1017  molecule  cm8  AT  300  K 

•  IN  Ar/H2  MIXTURES  [H2]  “lx  lO*8  molecule  cm"8  WILL 
QUENCH  BOTH  A r*  AND  H2(B) 


VLoz-bUb 


H2(B)  collision  frequency  as  a  function  of 

DISTANCE  FROM  THE  NOZZLE  THROAT 


250  m  DIAMETER  NOZZLE,  ?Hz  =  5  ATM 


1011 


10 10 


COLLISION 
FREQUENCY  109 

(s-1) 


io7L 

o.o 


4  0.8  1.2  1. 

EXPANSION  DISTANCE  (mm) 


A 

*  IV 

>1 


»VI 


Vi 

& 


Vub 7-^oc 


STABILIZATION  OF  H41 


LIFETIME  OF  H2(B)  -  H2  COMPLEX  10"10  -  10'12s 


H2(X)  +  hv  +  H2(B,v) 
H2(B)  +  H2  -►  H4* 


J’l  =  al 


H4  +  M  >  2H2  +  M 


H4  +  M  >  H4  +  M 


H4*  -►  H2(B)  +  H2 


H2( B)  H2  +  hv 


H4  LIFETIME  SHORT  -  STEADY-STATE  SOLUTION  APPROPRIATE 


d[H4]  ^  f  4>k2[H2]  k3[H2] 

~at ~  =  Ji[ri2j  1\2'[h'2t+-^*  ^3ih2:  +  k" : 


ESTIMATED  H4  FORMATION  EFF I C I ENC I ES/h2 ( B ) 
[H2]  =  lO1^  molecule  cm'2 


k 4  (s'1) 


Ha  EFFICIENCY 
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PROCEEDINGS  Or  THE  HIGH  ENERGY  DENSITY  WITTER  <HEDH> 
CONFERENCE  HELD  IN  R. .  <U>  RIR  FORCE  RSTRONRUTICS  LM 
EMMRDS  DFICII  HJ  LMJOERDRLE  ET  AL.  SEP  07 


UNCLASSIFIED  RFRL-CP-07-002 


DETECTION  OF  H4 


PHOTODISSOCIATION  TO  FORM  H2(B) 

X  <  500  nm 

DETECT  H2(B)  EMISSION  WITH  PMT 

FLASHLAMP  PHOTOLYSIS.  X  =  200  TO  800  nm 

LASER  PHOTOLYSIS.  X  *  337  nm.  266  nm 

CROSS  SECTION  UNCERTAIN:  10"18  TO  10‘20  cnfl2 

AT  10"20  cm2  CAN  DISSOCIATE  5  TO  10  PERCENT  WITH 
337  nm 

PHOTODISSOCIATION  TO  FORM  H3  +  +  H" 

X  <  670  nm 

COLLECT  TOTAL  IONS  FORMED  WITH  ELECTRODES 

POSSIBLE  CONSTRUCTION  OF  TIME-OF-FLIGHT  MASS 
SPECTROMETER 

IMPROVED  COLLECTION  EFFICIENCY  FOR  IONS  INCREASES 
EXPECTED  SIGNAL  BY  104 

PASSIVE  EMISSION 

QUASI -CONTINUOUS.  X  >  153  nm 

LOOK  WITH  PMT/VUV  MONOCHROMATOR 


WAVELENGTH 
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ABSTRACT 


Theoretical  Studies  of  Metastable  States  of 


Byron  H.  Lengsfield  III,  George  F.  Adams,  Cary  F.  Chabalowski, 
and  James  0.  Jensen 

US  ARMY  Ballistic  Research  Laboratory 
SLCBR-IB-I 

Aberdeen,  MD  21005-5066 


David  R.  Yarkony 
Department  of  Chemistry 
Johns  Hopkins  University 
Baltimore,  MD  21218 


We  have  recently  developed  a  series  of  theoretical 
methods  which  are  capable  of  calculating  the  quenching 
rates  of  metastable  molecular  states  due  to  nuclear 
and  relativistic  non-adiabaticity .  We  currently  have 
the  ability  to  analytically  compute  the  first  and 
second  order  nuclear  non-adiabatic .coupling  matrix 
elements,  <  V|#  n/JRI  Mf  >  and  l®7dfv  M*  >'  as  wel1  as 

to  compute  quenching  rates  arising  nrohi  "pin-orbit 
interactions.  The  ability  to  include  the  effect  of 
spin-spin  interactions  is  being  incorporated  in  our 
program  package  and  we  will  report  the  progress  attained 
in  that  effort. 

These  methods  have  been  applied  to  the  studies  of 
H4  and  He2,  and  we  will  report  the  results  of  that 
work.  For  H4  we  have  optimized  the  structures  of 
the  low-lying  states  at  the  multi-reference  Cl  level. 

The  lowest  state  in  the  pyramidal  (  C3V  )  configuration 
was  found  to  be  unstable  with  respect  symmetry  breaking 
nuclear  displacements.  The  first  excited  state  was  of  E 
symmetry  in  C3V  and  was  found  to  distort  to  a  Cs 
minimum  with  full  geometry  optimization.  Dipole  transition 
moments  and  non-adiabatic  coupling  matrix  elements  were 
calculated  at  the  minimum  of  the  excited  state. 

For  He2  we  have  considered  the  spin-forbidden 
radiative  transition  between  the  lowest  triplet  state 
and  the  ground  state,  — >  Y  •  We  have  also 

calculated  the  spin-orbitcoupling  matrix  elements  which 
are  needed  to  compute  the  quenching  of  the  lowest 
triplet  atomic  state. 
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Theoretical  Studies  of  Metastable 
States  of  Hei  and  H. 


Develop  Theoretical  Methods 
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Theoretical  Methods  Used  to 
Characterize  Potential  Energy 
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10,000  CSFs 


{Multi-Configuration  description  is  needed  for 
excited  states  and  for  systems  where  resonance 
is  important} 


See  p.  5792  of  J  Am  Chem  Soc,  108(19),  1986. 
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TABLE  II.  THEORETICAL  INTENSITIES  COMPUTED  WITH  HYBRID  BASIS  SETS 

(Intensity  in  KM/MOLE) 


* 


Spin-Orbit  interactions 


{Sun  over  states  is  avoided  by  solving  a  large  set  of  linear  equations 
(eqn  2.5a  of  reference).  This  means  large,  direct-CI  vavefunctions  can 
be  employed  and  the  slov  convergence  of  a  sum  over  states  perturbation 
expansion  is  avoided.  See  Yarkony,  J  Chem  Phys,  85,  7261(1986)) 


(Same  approach  is  being  developed  for  spin-spin  interaction) 


- NEXT  SLIDE - 

{We  can  exploit  work  on  derivatives  of  Cl  vavefunctions  to  compute 
non-adiabatic  matrix  elements) 

See  J  Chem  Phys,  81(10),  4549(1984). 

- NEXT  SLIDE - 

See  equation  2.8  of  J  Chem  Phys,  81(10),  4549(1984) 

(Two  terms  need  to  be  evaluated) 

See  equation  2.12a  and  2.12b  of  J  Chem  Phys,  81(10),  4549(1984) 


(Trace  of  a  one-particle  density  matrix  with  a  half-derivative  overlap 
integral) 


Non -Adiabatic  Coupling 

Born-Oppenheimer  Approximation 


vfSLm 


See  equation  2.13b  of  J  Chea  Phys,  81(10),  4549(1984) 


(Cl  gradient  evaluated  vith  a  transition  density  aatrix) 


-NEXT  SLIDE- 


See  Figure  1  of  Ginter,  H.L.  and  Battino,  R.,  4470. 


(A  large  number  of  curve  crossings  must  be  considered.  Relativistic 
interactions  can  result  in  avoided  crossings.) 
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Theoretical  Studies  of  He 

Input  Parameters  to  Liquid  Simulation 


Lifetime  of  the  *2^  State  of  He 
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Crystal  City,  VA 

The  Influence  of  Condensed  Phase  on 
Metastability 

by 

P.  K.  Swaminathan*,  C.  S.  Mur  thy,  B.  C.  Garrett,  and  M.  J.  Redmon 
Chemical  Dynamics  Corporation 
9360  Pennsylvania  Avenue  #106 
Upper  Marlboro,  MD  20772 

Condensed  phase  influence  on  HEDM  candidates  plays  a  critical  role  in 
determining  the  lifetimes  and  storage  conditions,  thereby  ultimately 
influencing  their  practical  usefulness.  We  have  developed  a  comprehensive 
theoretical  methodology  involving  quantum  chemistry  input,  modern 
computer  simulation  techniques  (such  as  Monte  Carlo,  Molecular  Dynamics, 
and  Generalized  Langevin  Dynamics),  and  semidassical  eikonal  description  of 
electronic  inelasticity  to  investigate  various  aspects  of  the  underlying 
microscopic  energy  transfer  mechanisms. 

We  will  show  illustrative  results  obtained  with  the  new  methodology 
for  the  well  known  examples  of  liquid  Helium  metastables,  so  far  employed 
in  this  research  as  convenient  prototypes.  The  experience  gained  in  modeling 
metastability  in  these  cases  has  laid  the  foundation  for  future  studies  on 
other  potential  HEDM  candidates. 
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Our  Role: 

Develop  Methods  for  Elucidation  of 
Quenching  Mechanisms  in  Gas 
and 

Condensed  Phases 


?  HEDM  ? 

Scenario: 


HEDM  s  involve  energetic  sfecies,  usually 
in  an  ELECTRONICALLY  EXCITED  state 

Therefore,  the  chemistry  and  quenching 
mechanisms  involve 

ELECTRONICALLY  NONADIABATIC 
Collision  Dynamics 

Condensed  Phase  materials 
low  molecular  weight  inert  species 


Temporary  Prototypes: 

Helium  Metastables  He  (3S  )  an(* 

He*  (a  -  state)  in  Helium  Liquid 
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Technical  Approaches 

Gas-Phase 

Self-Consistent  Eikonal  Method  (SCEM) 

Electronic  amplitudes  coupled  to 
Nuclear  Trajectories  via  an  Ehrenfest 
Effective  Potential 

A  Self-Consistent  Semiclassical 
Description  practical  for  polyatomics 

Condensed  Phase 

Computer  Experiments:  Monte  Carlo, 

Molecular  Dynamics 

Few-  body  Models  • 

SCEM  +  Generalized  Langevin  Equation  (GLE) 
based  STOCHASTIC  MODELS 

System  +  Heatbath  Decomposition  with 
one-time  heatbath  parameterization  via 
computer  experiments 


He-He  Effective  Pair  Potentials 


r’(r’=r/ a) 


U(r’)/k  (K) 


Liquid  Structures  of  Model  Helium 
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<Y(o)-Y(t)>/<v(o)2> 


e  I  i  u  m 


14.3  GPa 
12.7  GPa 
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0.3  GPa 


<V(0)«V(t)>/<V(0)2> 


U(r)/hartree 
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0  04 


0  01 


-0  01 


-0  03 


-0  06 


-0  08 
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OPTIMIZED  TRIAL  FUNCTIONS  FOR  QUANTUM 
MONTE  CARLO  STUDY  OF  H4* 

Sheng-yu  Huang, 1  Zhiwei  Sun,^  and  William  A.  Lester,  Jr.^ 

Materials  and  Chemical  Sciences  Research  Division 
Lawrence  Berkeley  Laboratory 
University  of  California 
Berkeley,  CA  94720 

In  recent  years  the  fixed-node  quantum  Monte  Carlo  method  (FNQMC) 
has  been  successfully  applied  to  the  calculation  of  electronic  energies  of  small 
molecules.  The  attractive  feature  of  this  method  is  that  with  a  modest  basis  set 
and  often  a  single-Slater  determinant  one  typically  recovers  over  90%  of  the 
correlation  energy. 

Because  the  simulation  is  dominated  by  a  trial  function  that  fixes  the 
nodal  positions,  a  good  trial  function  is  essential  in  order  to  obtain  an  accurate 
result.  It  is  often  chosen  to  be  an  ab  initio  approximate  Hartree-Fock  wavefunc- 
tion  in  which  case  the  many-electron  nodes  are  established  without  consideration 
of  electron  correlation.  Introduction  of  an  explicit  correlation  factor  in  the  Jas- 
trow  form  (exponential)  improves  the  variational  energy,  but  this  nodeless  func¬ 
tion  does  not  improve  the  nodes  of  the  system.  Ideally,  the  Slater  determinant 
needs  be  re-optimized  in  the  presence  of  the  correlation  factor. 

We  have  developed  an  algorithm  to  optimize  trial  function  parameters 
(linear  coefficients,  non-linear  exponential  factors,  and  parameters  of  the  correla¬ 
tion  function)  based  on  a  random  walk  procedure.  In  order  to  obtain  accurate 
energies  for  excited  states,  we  introduce  a  projection  operator  method  based  on 
group  theory  to  constrain  the  optimization  process.  It  effectively  distinguishes  all 
the  states  of  a  given  symmetry.  We  also  apply  the  Young  tableaux  of  permuta¬ 
tion  group  theory  to  facilitate  the  treatment  of  fermion  properties  and  muitiplcts. 

Using  an  optimized  trial  function,  we  are  studying  pyramidal  H4  at  the 
maximum  ionicity  excited  state1  (MIES)  and  other  geometries  relevant  to  under¬ 
standing  decomposition  pathways  in  this  system.  The  calculated  potential  ener¬ 
gies  show  0.6-0.8  eV  lowering  compared  to  a  previous  study  .  In  a  separate 
multi-configuration  ab  initio  study  we  have  determined  that  there  is  a  pathway 
without  a  barrier  for  H2(B)  approach  to  ground  state  H2(X).  Work  is  in  progress 
to  fully  characterize  the  first-excited  state  of  the  H4  system  and  its  coupling  to 
the  ground  state,  and  will  be  discussed. 

*  This  work  was  supported  by  the  U.S.  Air  Force  Rocket  Propulsion  Laboratory  (AFRPL) 
through  agreement  with  the  Department  of  Fnergy  under  Contract  No.  DK-AC03-76SF00098. 

1  Also,  Department  of  Chemistry,  University  of  California,  IJerkelcy,  CA  9-1720 

*  Permanent  address:  Institute  of  Mechanics,  Academia  Sinica,  People’s  Republic  of  China 

1  C.  A.  Nicolaides,  G.  Theodorakopoulos,  and  I.  D.  Petsalakis,  ,1.  Chein.  Phys.  80  (4),  1705 
(1984).  ' 


QMC  AND  MCSCF  STUDY  OF  TETRAHYDROGEN 

j.  f 

Sheng-Yu  Huang,'  Zhiwei  Sun,  and  William  A.  Lester,  Jr. 

Materials  and  Chemical  Sciences  Division 
Lawrence  Berkeley  Laboratory 
University  of  California,  Berkeley 
Berkeley,  California  94720 

I.  MCSCF  Pilot  Study  of  Trigonal  Pyramidal  Ground-  and  Excited-State  H^. 

We  have  demonstrated  that  MCSCF  (multiconf iguration  self-consistent  field) 
wave  functions  as  trial  functions  for  FNQMC  (fixed-node  quantum  Monte  Carlo) 
calculations  can  recover  -10011  of  the  correlation  energy  of  ground  states  and 
~95%  of  the  correlation  energy  of  excited  states.^  Using  an  extended  (triple 
zeta-plus-polarization  (TZP))  basis  set,  MCSCF  calculations  were  carried  out 
for  the  MIES  (maximum  ionicity  excited  state),  see  Fig.  1.  They  confirmed 
the  characteristics  depicted  in  Fig.  2  for  the  C3y  pyramidal  structure 
including  the  minimum  energy  geometry,  the  doubly  degenerate  E  ground  state  as 
H  separates  from  for  R  >  3.8  a.u.,  and  the  nondegenerate  A  ground  state 
for  R  <  3.8  a.u.,  consistent  with  an  E-  to  A-state  crossing  at  R  =  3.8  a.u. 

For  C$  symmetry  the  avoided  crossing  was  obtained.  Figure  3  displays  the 
results  of  the  MCSCF  pilot  computations.  These  curves  lie  above  those  of 
Nicolaides,  Theodorakopoulos,  and  Petsalakis  (NTP)  and  reflect  the  better 
capability  of  the  ab  initio  MRD-CI  wave  functions  compared  to  the  compact  MCSCF 
trial  functions  needed  here  as  FNQMC  trial  functions.  The  key  comparison, 
however,  will  be  with  the  FNQMC  results  obtained  with  the  MCSCF  functions.  We 
discuss  this  point  in  Sec.  III. 

II.  MCSCF  Pilot  Study  of  the  H?(X)  ♦  H?(B)  *  H4  MIES  Pathway. 

Calculations  were  carried  out  for  H^B)  approach  to  H^(X)  in  which 

these  molecules  are  contained  initially  in  perpendicular  Dianes  that  bisect 


each  other.  Figure  4  displays  this  arrangement  and  summarizes,  in  the 
accompanying  table  a  sequence  of  geometry  changes,  labeled  I -VI I,  that  carry 
the  system  from  the  asymptotic  region  to  the  neighborhood  of  the  MIES 
configuration.  The  table  also  contains  the  energy  lowering  associated  with 
the  steps  that  are  plotted  in  Fig.  5.  It  is  noteworthy  that  no  energy  barrier 
is  encountered  along  this  path.  Further,  Fig.  6  shows  that  when  H(B)  is 
displaced  from  bisecting  H^X)  in  step  II,  charge  transfer  immediately 
occurs.  (Similar  behavior  has  been  found  for  the  related  H^B)  +  He  system 
in  an  independent  study  by  WAL.) 

III.  FNQMC  Study  of  the  Ground-State  Pyramidal  Structure  for  C^v  Symmetry 
FNQMC  calculations  using  the  MCSCF  trial  functions  discussed  in  Sec.  I 

yield  ~  0. 6-1.0  eV  energy  lowering  compared  to  the  results  of  NTP  and  are 
presented  in  Fig.  7.  Such  a  large  change  was  not  anticipated  and  so  it  was 
important  to  test  the  validity  of  this  finding.  To  this  end  a  configuration 
interaction  calculation  including  all  single  and  double  excitations  (SDCI) 
using  the  MCSCF  pilot  study  basis  set  was  carried  out  at  R  =  3.4  a.u.  The 
energy  was  0.32  eV  lower  than  NTP's  value  and  is  consistent  with  the 
improvement  expected  based  on  studies  of  other  systems.  The  FNQMC  results  of 
Fig.  7  obtained  using  a  new  trial  function  optimization  algorithm,  see  Sec. 

IV,  are  generally  an  improvement  over  those  of  Fig.  7  obtained  using  MCSCF 
trial  functions. 

IV.  FNQMC  Study  of  the  Excited  State 

These  calculations  provide  the  severest  test  of  the  FNQMC  approach  because 


of  the  lack  of  knowledge  of  the  accuracy  of  the  excited  state  trial  function 
needed  to  provide  a  nodal  description  that  assures  orthogonality  to  the  ground 
state  of  the  same  (A1)  symmetry  (Cs).  MCSCF  calculations  close  to  the 
avoided  crossing  suffered  from  root  flipping.  Despite  the  use  of  familiar 
MCSCF  strategies  to  address  the  problem,  it  could  not  be  resolved. 

The  MCSCF  convergence  problem  had  been  encountered  earlier  for  He  +  H^B) 
but  resolved  by  the  use  of  an  ab  initio  Cl  method.  Interest  here  in  using  the 

3 

FNQMC  method  led  us  to  develop  a  method  for  trial  function  construction 
that  avoids  the  MCSCF  procedure.  The  approach  amounts  to  the  introduction  of 
parameter  optimization  in  the  random  walk  process.  Using  group  theory  a 
projection  operator  is  constructed  and  used  to  constrain  the  wave  function  to 
have  the  symmetry  properties  of  the  state  of  interest.  This  method  has  been 
employed  in  the  present  effort  to  generate  the  A-state  FNQMC  results,  denoted 
A(FNQMC) ,  of  Fig.  7.  Calculations  using  this  method  are  in  progress  to 
complete  the  E-state  curve,  labeled  E(FNQMC),  in  Fig.  7  for  3.8  <  R  <  6.0  a.u. 

V.  SA-MCSCF  Study  of  Distorted  Geometries 

Calculations  have  been  carried  out  using  the  state  averaged  (SA)-MCSCF 
method  to  develop  trial  functions  for  a  QMC  study  of  the  topography  of  the 
ground-  and  excited-state  potential  energy  surface  (pes)  in  the  region  of  the 
symmetric  geometry  of  the  MIES  determined  by  NTP.  Figure  8  presents  the 
coordinate  system.  Fig.  9  provides  perspective  views,  and  Fig.  10  (R  =  4.0 
a.u.)  ana  Fig.  11  (R  =  3.8  a.u.)  show  contour  maps  of  the  pes  in  the  MIES 
region.  Figures  10  and  11  show  that  the  ground  state  has  a  saddle  point  at 
smaller  R  than  the  minimum  of  the  excited  state  and  that  both  features 


correspond  to  an  isosceles  triangle  base  for  the  MIES  system.  Further 
geometry  optimization  is  explored  in  Fig.  12  which  plots  potential  energy  as  a 
function  of  d,  the  displacement  of  H  towards  the  base  of  the  H,  isosceles 
arrangement.  The  minimum  for  the  excited  state  is  found  for  d  =  0.1  a.u. 

VI.  Nonadiabatic  Coupling 

The  stability  of  the  MIES  system  is  dependent  on  nonadiabatic  coupling 
(NAC)  to  the  ground  state.  Although  the  focus  of  this  study  is  characteriza¬ 
tion  of  the  region  of  pes  of  the  MIES,  tne  need  to  ascertain  first  the 
stability  of  the  excited  system  is  a  high  priority  here  because  of  the  effort 
associated  with  the  trial  function  optimization  method  and  the  expense  of  QMC 
calculations  with  small  statistical  variances.  To  calculate  the  NAC  matrix 
elements  we  took  advantage  of  the  simplifications  made  possible  by  the  use  of 
normal  mode  distortions  (Fig.  13)  from  the  highly  symmetric  C3y  symnetry. 

Our  calculations  show  that  nuclear  displacements  associated  with  modes 
and  lead  to  strong  coupling  between  the  A  and  E  states.  Relatively 
weaker  coupling  is  found  in  and  Q^.  The  two  totally  symmetric  modes 
(Q3  and  Qg)  give  no  contribution. 

The  size  of  the  NAC  matrix  elements  connecting  tne  A'  state  to  the  E-state 
components  (C3v  notation  is  used  here  to  indicate  parentage)  are  tabulated 
in  Fig.  14  for  the  minimum  energy  geometry  of  the  A  state.  These  results 
establish  that  there  is  strong  coupling  between  these  states.  Further 
computational  study  is  needed  to  confirm  these  predictions  and  estimate 
reliabily,  the  excited  state  lifetime. 
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Figure  Captions 


Figure  1.  MIES  coordinate  system  (assuming  C^v  symmetry):  R  is  the  distance 
from  Ha  to  plane  of  H^H  equilateral  triangle  of  side  r. 

Figure  2.  Ground-  and  excited-state  potenial  energy  curves  for  trigonal 

pyramidal  geometry:  upper  panel  (C^v  symmetry),  lower  panel  (C^  symmetry). 

Figure  3.  Potential  energy  curves  from  MCSCF  trial  functions  for  trigonal 
pyramidal  geometry  and  C^v  symmetry  (r  =  1.70  a.u.).  State  designations 
for  symmetry  are  indicated  in  parentheses.  The  lack  of  coincidence 
of  the  degenerate  E  curves  reflects  their  calculation  in  lower  symmetry 
and  provides  an  indication  of  MCSCF  convergence. 

Figure  4.  Geometries  used  to  calculate  an  approach  of  H^B)  to  ( X ) . 

Figure  5.  Plot  of  total  energies  corresponding  to  the  pathway  presented  in 
Figure  4. 

Figure  6.  Change  of  charge  distribution  between  H  and  H,  as  a  function 

d  D 

of  d,  the  displacement  of  the  midpoint  of  H?(B)  above  the  plane  of  H?(X). 

Figure  7.  Potential  energy  of  trigonal  pyramidal  versus  R.  The  A  (varia¬ 
tional)  results  were  determined  using  a  new  optimization  method  for  QMC 
trial  functions. 

Figure  8.  MIES  coordinate  system  used  for  geometry  variation,  d  is  the 
magnitude  of  the  displacement  o*  Ha  along  the  bisector  of  the  HcHd 
side  of  length  b. 

Figure  9.  Potential  energy  surfaces  for  the  ground-  and  excited-state  sepa¬ 
rately  for  easy  visualization  and  in  computed  relative  positioning.  Note 
the  minimum  of  the  excited  state  and  saddle  point  of  the  ground  state. 


Figure  10.  Potential  energy  contour  map  for  R  =  4.0  a.u.  as  a  function  of  a 


and  b;  see  Fig.  8. 


Figure  11.  Same  as  Fig.  10  for  R  =  3.8  a.u. 


Figure  12.  Potential  energy  curves  for  selected  R  (in  a.u.)  vs  d;  see  Fig.  8. 


Figure  13.  Normal  mode  displacments  of  in  C^v  symmetry. 

Figure  14.  Cl  contribution  to  coupling  matrix  elements  for  normal  mode 


displacement,  see  Fig.  13. 


ENERGY  (a. 


Geometries  and  energies  for  the  approach  of  H2(B)  to  H2(X) 


step 

R; 

Lj  U 

L2(+) 

Lo(-) 

D 

energy  lowing 

I 

5-2.2 

1.40  2.43 

1.215 

1.215  0.0 

-0.01 

n 

2.2 

1.40  2.43 

1.215-1.03 

1.215-0.5 

0.715 

-0.20 

m 

2.2 

1.60  2.63 

1.03-2.13 

0.5 

0.865 

-0.20 

IV 

2.2-1.75 

1.60  2.63 

2.14 

0.5 

0.865 

-0.20 

V 

1.75-1.55 

1.60  2.73 

2.23 

0.5 

0.015 

-0.46 

VI 

1.55 

1.60  3.13 

2.23-2.63 

0.5 

1.115 

-0.60 

vn 

1.50 

1.70  3.23 

2.63-2.73 

0.5 

1.165 

-0.11 

TOTAL  ENERGY  LOWERING  FOR  THE  SEVEN  STEPS  IS  5.75 
ev. 


*  Distances  in  a.u.;  energies  in  ev. 

R'  -  distance  between  the  midpoint  of  H2(X)  and  the  point  where  H^B)  meets 
the  H2(X)  plane. 

L|  -  length  of  H2(X). 

4  -  length  of  H^B). 

Li(+)  -  length  of  H2(B)  above  H^X)  plane. 

Li(-)  -  length  of  H^B)  below  H^X)  plane. 

D  -  magnitude  of  shift  of  H2(B)  midpoint  (  above  (+)/beiow(-) )  plane  of  H2(X). 
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SHIFT  OF  MIDPOINT  OF  ll,(D)  FROM  BISECTOR  PLANE 


POTENTIAL  ENERGY  OF  PYRAMI  DAL  H 


DISTANCE  BETWEEN  H  AND  H,  fa.u.I 


Fig.  1  excited  state 
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CONTOUR  MAP  (R— 4.0  a.u.) 
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PI 

Coupling  matrix  elements  D  by 


a  finite  difference  method  (  extrapolated  from 
<5Q= 0.002,  0.001,  and  0.0005  ). _ 


<E,  |  A'  > 

<E,  |  A'  > 

Qx 

-2.0112(64) 

0.0015(0) 

q2 

0.0017(0) 

2.0127(98) 

Qa 

0.0000 

0.0000 

Q< 

-0.3904(389) 

-0.0010(3) 

Qs 

-0.0039(13) 

-0.3113(196) 

Qa 

0.0000 

0.0000 

1.  Ej  and  E2  are  doubly  degenerate  states  where  El  is 
symmetric  and  E2  is  antisymmetric. 


Abstract  for  HEDM  Conference 
Spectroscopy  of  Polyatomic  Hydrogen  Ions 

Takeshi  Oka 

Infrared  laser  spectroscopy  of  the  v^-fundamental  band  of  H3+  and  the 
extension  of  work  towards  higher  energy  state  and  isotopic  species  will  be 
reported.  The  chemical  and  physical  dynamics  of  this  ion  in  gaseous  discharge 
will  be  discussed.  Our  plan  to  observe  spectra  of  and  its  cluster  ions  in 
condensed  phase  plasma  will  be  discussed. 
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MFASURFMFNT  Of  rill  (  I IARGE  rRANSFFR  RA  IT  CONSTANT  COR.  D34  ,  Sil  !/, 

P.  D.  I  iaaland  and  A.  Garscadden 
Air  Force  Wright  Aeronautical  Laboratories 
Wright-Patterson  Air  Force  Base 


Studies  of  dissociative  charge  transfer  in  silane  collisions  with  SiFl2T  and  Sih3  + 
ions  have  indicated  that  I  t-  is  transferred.  In  silane  deposition  reactor  plasmas,  H,+ 
and  H2+  also  are  formed  and  are  expected  to  rapidly  form  I  13  +  .  The  reaction  of  I  13 + 
on  silane  is  therefore  of  interest.  In  order  to  resolve  ambiguities  in  the  reaction,  the 
measurements  were  made  using  03+.  The  experiment  utilised  a  modified  Nicolet 
Fourier  Transform  Mass  Spectrometer.  A  mixture  of  deuterium  and  silane  was  ionized 
by  an  electron  beam.  The  silane  ions  were  then  ejected  from  the  trap.  The  deuterium 
molecular  ions  react  rapidly  with  background  deuterium  to  form  D3R  The  F>3+  reacts 
with  the  background  silane  to  give  SiH3+  and  neutral  products.  This  rate  constant  is 
estimated  as  approximately  2x  lO'^cm^s-  I . 


MEASUREMENT  OF  THE  CHARGE  TRANSFER 
RATE  CONSTANT  FOR  BX  +  SiH* 
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ABLY  ASSISTED  BY  MR.  JIM  BARNHART 


H'  TRANSFER  OBSERVED  IN  DOUBLE  RESONANCE 
EXPERIMENTS 


HENIS,  STEWART,  TRIPOD!  AND  GASPAR,  J.  CHEM.  PHYS.  57,  389  (1972) 


FOURIER  TRANSFORM  MASS  SPECTROMETRY 


6.  EXPERIMENTAL  CONFIGURA 

laser  PLASMAS,  electro* 


DEUTERIUM  TRIMER  ION 


MASS  IN  A.M.U. 


DEUTERIUM  TRIMER  ION 


MASS  IN  A.M.U. 
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EXPLOSIVE  DECOMPOSITION  OF  FLUORINE  AZIDE  FILMS 


D.J.  Benard 


Rockwell  Science  Center 
1049  Camino  dos  Rios 
Thousand  Oaks,  CA  91360 


ABSTRACT 


The  molecule  FN^,  known  as  fluorine  azide,  can  be  viewed  as  the 


metastable  species  NF(a  A)  bound  with  ^ .  The  nature  of  the  binding 


will  be  discussed  in  terms  of  the  potential  surfaces  and  their 


correlations  to  the  various  excited  electronic  states  of  Nj  and  NF. 


Experimental  data  will  also  be  presented  in  the  form  of  emission  spectra 


obtained  from  the  laser  initiated  detonation  of  thin  FN^  films.  These 


data  are  expected  to  shed  some  new  light  on  methods  to  stabilize  high 


concentrations  of  metastable  molecules  which  may  be  useful  as  high  impulse 


rocket  propellants.  The  use  of  FN^  as  a  starting  material  for  the 


generation  of  high  concentrations  of  Nt(a  A),  to  power  short  wavelength 


laser  systems,  will  also  be  discussed. 
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DISTORTED  MOLECULE  CALCULATION 
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KEY  ELEMENTS  OF  SWCL 
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ALTERNATIVE  GENERATION  SCHEME 


FLOWTUBE  EXPERIMENT 


REACTION  STEPS 
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Electronic  Structure  and  Stability  of  Energetic  Chemical  Species 

H.  H.  Michels  and  J.  A.  Montgomery,  Jr. 

United  Technologies  Research  Center 
East  Hartford,  CT  06108 

ABSTRACT 

Theoretical  quantum  mechanical  calculations  have  been  carried  out  for 
several  potentially  high  energy  chemical  species.  Our  research  objectives 
are  to  identify  those  light  element  molecular  species  that  have  sufficient 
internal  energy  to  be  useful  in  advanced  chemical  propulsion  systems.  A 
parametric  study  of  the  deliverable  specific  impulse  as  a  function  of 
available  reaction  energy  and  propellant  weight  indicates  that  primary 
consideration  should  be  given  to  molecular  structures  that  can  be  formed  from 
hydrogen  through  boron  and  that  have  a  molecular  weight  of  40  or  less. 

Calculations  to  date  have  been  performed  on  three  classes  of  light 
element  compounds:  C^v  structures  such  as  and  Li^H,  azide  structures 
such  as  N^F  and  0^02  and  cyclic  boron  structures  such  as  B^H^  and 
B2H2NH.  We  find  that  the  ground  state  potential  energy  minimum  found  for 
in  C^v  symmetry  corresponds  to  a  saddle  region  rather  than  a  stable 
bound  state  and  that  distortion  via  vibrational  or  rotational  modes  leads 
monotonically  to  dissociation  into  two  H2  molecules.  In  contrast,  Li^H 
appears  to  be  chemically  stable  in  C^v  symmetry.  Our  studies  of  asymmetric 
dinitrogen  dioxide  (  indicate  a  stable  azide-like  structure  in 

Cs  symmetry  for  the  ground  *A'  state,  with  an  indicated  heat  of  formation 
of  +430  kjoule/mol.  This  molecule  is  a  very  attractive  candidate  as  an 
advanced  oxidizer.  Our  preliminary  calculations  for  boron  compounds  indicate 
that  B^Hj  is  unstable  in  symmetry  but  that  iminodiborane 
(B2H2NH)  is  stable  in  C2V  symmetry.  The  electronic  structure, 
vibrational  analysis  and  thermodynamic  stability  of  these  compounds  will  be 
described . 


SPECIFIC  IMPULSE,  Isp  (sec) 


luviii'vm’uw  wi 


wwrvTi?iviw  v*  \r*  y-*  v»  «  ■ 


SPECIFIC  IMPULSE  OBTAINABLE  FROM  FUEL/OXIDANT  ENERGY  RELEASE  WITH 

HYOROGEN  AS  A  WORKING  FLUID 

HYDROGEN  CONTENT  varied  to  OfMiMtfE  IMPULSE  (1000  -><  ’  ps^i 


•  6  «  3(. 


CATEGORIES  OF  ENERGETIC  SPECIES 


1)  Ordered  (or  metallic)  forms  of  matter  which  at  conventional 
temperatures  and  pressures  are  molecular  solids.  Metallic 
hydrogen  and  ammonium  are  examples 

2)  Free  radicals,  either  atomic  or  molecular  in  nature,  that 
are  stabilized  against  recombination  by  condensation  and 
immobilization  on  noble  gas  matrices  held  at  very  low  (4  K) 
temperatures.  Atomic  hydrogen,  beryllium  monohydride  (BeH) 
and  the  hydroxyl  (OH)  radical  are  examples  in  this  category. 

3)  Electronically  excited  metastable  atoms  or  molecules  that 
are  stabilized  against  radiative  decay  by  externally  applied 
electric,  magnetic  or  radiation  fields,  or  by  an  intrinsic 
field  within  a  condensed  phase  of  such  material  which  acts 
as  a  barrier  to  their  decay.  Metastable  atoms,  such  as 

o  p 

He*(  S)  or  N*(  D),  and  metastable  molecules,  such 

as  He  *(a3X  +  or  5£  +),  HeH  *  ^A  II)  and 
2-j  .  u  g  2 

H  ‘  (B  £u),  are  examples  in  this  category  that  would 
constitute  significant  energy  sources  (>  1000  kjoule/mo!). 

Other  energetic  metastable  species  such  as  O  (a  dg),  NF 
(a Id)  and  Mg*(3P),  with  long  radiative  lifetimes,  may 
also  be  useful  as  temporary  energy  storage  media. 

4)  Strained  molecules  with  high  positive  heats  of  formation. 

This  is  an  important  class  of  molecules  which  are  really 
isomeric  forms  of  conventional  chemical  species  but  are 
structurally  prevented  from  relaxation  to  their  lowest 
energy  stale  by  significant  activation  barriers  on  their 
potential  energy  surface.  Examples  in  this  category  include 
cycloproplyene  (C^H^),  aziridine  (C^^NH),  azetidine  (C^H^NH), 
hydrogen  azide  (HN  )  and  fluorine  azide  (FN  ). 
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A.  HIGH-ENERGY  LIGHT  ELEMENT 
MOLECULAR  STRUCTURES 


c)  LijH  (C3v)  d)  Li4  (C3v) 


H4(C3v)  calculations 
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LUH  CALCULATIONS 


Energies  (hartrees),  distances  (A  ),  frequencies  (cm’ 
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B.  AZIDE-TYPE  STRUCTURES 


BM13OTX.0OH 


ENERGETICS  OF  THE  CHEMICAL  SYNTHESIS 


fn3  dissociation  pathway 


FN-N2  reaction  coordinate  (A) 


FLUORINE  AZIDE  (FN3) 
FREQUENCY  ANALYSIS 
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OPTIMIZED  GEOMETRIES 
OF  AZIDES 


Experimental  bond  length  (A) 


1.098 

1.128 


H-N 

0-0 


1.041 

1.216 


1.308 

1.151 


1.281 

1.262 


Compound 

Theory 

Rl  (A) 

R2  (A) 

R3  (A) 

$ 

a  (deg)  [3  (deg)  >; 

hn3 

SCF 

1.0055 

1.2381 

1.0987 

108.181 

173.815  | 

fn3 

scf 

1.3820 

1.2536 

1 .0995 

104.315 

174.108  V, 

MP2 

1.4309 

1.2799 

1.1521 

103.765 

171.803 

FNCO 

SCF 

1.3737 

1.2387 

1.1354 

109.846 

173.235  3 

MP2 

1.4185 

1.2622 

1.1765 

110.717 

168.914  % 

a-N202 

SCF 

1.7574 

1.2024 

1.0844 

103.966 

179.506  £ 

179.488  i 

179.360  * 

MP2 

1.5305 

1.2272 

1.1548 

103.591 

CISD 

1.5815 

1.2238 

1.1072 

102.840 

FNBF 

SCF 

1.2993 

1.2068 

1.2871 

180.000 

180.000  3 

MP2 

1.3394 

1.2491 

1.3069 

0 

154.650 

169.291  $ 
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c.  HIGH  ENERGY  BORON 
STRUCTURES 

H 

I 


B 


a)  B3H3  (D3h)  b)  B2H20  (C2v) 

cyclotriborane  1,  2  epoxydiborane/diborylene  oxide 

c:  ethylene  oxide 


„/V. 

\/ 

B 

I 

H 

d)  B3H3NH  (C2v) 
triboryleneimine 
c:  azetidine/trimethyleneimine 


H 

I 


N 


c)  B2H2NH  (C2v) 
iminodiborane 
c:  aziridine/ethyieneimine 


CYCLIC  BORON  COMPOUNDS 

Cyclotriborane  (B3H3) 


Rgg  =  1 .626  A  Rgg  =  1 .734  A 

Rgl^  =  1 .232  A  Rgl^  =  1.178A 

E  =  -  75.038853,  stable  E=  -  75.174458,  unstable  (2) 


Iminodiborane  (B2H2NH) 

H 

I 


N 


rBH  =  1.174A  L  BNH  =  145.3° 

RNH  =  0.989  A  E=  -  104.981636 

RBN  =  1.423  A 


B1 _ A1 _ ^2 _ *2 

700  818  921  923 


Frequencies  (cm~  "*) 


989  1000  1110  1268 


A1  B2  A1  A1 
1423  2678  2916  3965 
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CONCLUSIONS 


1.  Tetrahydrogen  (H4),  LiH3  and  Li4  are  not  stable  molecules  in  their 
ground  state  in  C3v  symmetry. 

2.  Li3H  is  stable  in  its  ground  state  in  C3v  symmetry.  Additional 
calculations  of  the  energetics  of  Li3H  are  in  progress. 

3.  All  of  the  azide-like  structures:  FN3,  FNCO,  a-N202,  and 
FNBF,  are  vibrationally  stable  in  the  lowest  1A'  state  with 
predicted  high  positive  heats  of  formation.  FN3  also  has 
a  stable  low-lying  3A  state. 

4.  B3H3  has  no  stability  as  a  D3h  structure.  The  anion,  B3H3~, 
is  stable  as  predicted  by  Lipscomb's  rules.  B2H2NH  is  stable 
as  a  C2v  structure  with  symmetry  ^.j.  Preliminary  thermo¬ 
dynamics  indicate,  however,  that  B2H2NH  is  not  very  energetic. 

5.  Best  systems  studied  to  date:  a-N202  and  FNBF 

6.  Areas  for  futher  study  Hp.  LimHn  structures,  azides  and 
further  simple  boron  compounds 


METASTABLE  MOLECULAR  FUELS: 
THEORETICAL  STUDY  OF  ION-PAIR  STATES 
LOW-LYING  SURFACES  OF  H3O 

Roberta  P.  Saxon 
Dahbia  Talbi 


SRI  International 


OBJECTIVE 


Predict  new  energetic  metastable  molecular 
species  that  do  not  decay  by  radiation, 
tunneling,  or  other  means  when  isolated  in 
vacuum. 


PROPOSAL 

Ion-pair  bonding  between  stable  negative  ions 
may  lead  to  high-energy  (locally)  bound  states 
that  have  not  been  studied  previously. 


Example:  H4 


METASTABLE  MOLECULAR  FUELS: 


THEORETICAL  STUDY  OF  ION- PAIR  STATES  -  LOW-LYING  SURFACES  OF  H-jO 


Roberta  P.  Saxon  and  Dahbia  Talbi 
SRI  International 


Menlo  Park,  California  9402b 


Metastable  molecular  fuels,  long-lived  molecular  species  that  do 
not  decay  by  radiation,  tunneling,  or  other  means,  when  isolated  in 
vacuum,  have  been  proposed  as  the  basis  for  possible  new  propulsion 
schemes.  The  recent  prediction  of  an  energetic  excimer  state  of 
suggests  the  possibility  of  a  whole  series  of  molecules  bound  by  the 
Coulomb  attraction  between  a  stable  positive  ion  and  a  stable  negative 
ion.  The  fate  of  an  ion-pair  species,  once  formed,  will  depend 
sensitively  on  the  details  of  the  potential  surfaces  that  govern  decay 
processes  such  as  optical  transitions,  predissociation,  and  internal 
conversion.  This  theoretical  program  is  devoted  to  examining  ion-pair 
species  composed  of  first-row  atoms  that  could  serve  as  candidate 
fuels.  In  this  paper,  CASSCF-FOCI  calculations  on  the  low- lying 
doublet  and  quartet  potential  surfaces  of  the  H-jO  system  are  reported. 

Consistent  with  our  model  of  an  ion-pair  state  as  an 
equilateral  triangle  with  an  0*  located  above  the  center  of  the 
triangle,  our  calculations  to  date  have  concentrated  on  C^v 
geometries.  A  DZP  basis  set  augmented  by  diffuse  s  and  p  functions 
on  0  was  used.  The  correlation  diagram  linking  the  H^O  states  studied 
here  to  states  of  the  separated  fragments  has  been  established.  The 
connection  between  this  work  and  the  portions  of  the  ground  state 
(doublet)  potential  surface  considered  by  previous  workers^ i.e.  the 
H^O  local  minimum  and  the  H2  +  OH  -*  H2O  +  H  transition  state  has  also 
been  explored.  While  portions  of  the  lowest  potential  surface 
corresponding  to  the  ion-pair  H^-O’  can  be  identified,  at  the  minimum, 
energy,  the  lowest  state  may  be  described  as  H^0+  with  an  electron  in 
an  oxygen  Rydberg  orbital,  in  agreement  with  the  previous  report". 

This  conformation  is  known  to  dissociate  to  H2O  +  H.  There  is  no 
barrier  inhibiting  conversion  of  the  ion-pair  geometry  to  the  Rydberg 
geometry.  Within  the  restricted  C^  geometry,  however,  a  stable  ion- 
pair  state  with  H-H  distance  very  similar  to  that  in  the  ion  is 
observed  on  the  first  excited  doublet  potential  surface.  Investigation 
of  this  surface  in  unrestricted  geometries  is  underway. 


^Vork  supported  by  Air  Force  Flight  Test  Center  under  contract 
F04611-86-C-0070 
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S.  P.  Walch  and  T.  H.  Dunning,  J  Chem  Phys.  22.  -  >03  1980). 

K.  S.  E.  Niblaeus,  B.  0.  Roos ,  and  P.  E.  M.  Sit  gbahr.  Chem.  Phys. 
21,  207  (1977) . 


SELECTION  OF  SPECIES 


Maximum  specific  impulse 


heat  of  reaction 
mass  of  products 


low  molecular  weight 

H3+  positive  ion  -  1st  row  negative  ion 


optimal  electron  affinity  -  a  tradeoff 

large  EA  =>  stable  negative  ion 
small  EA  =»  high-energy  minimum 


ELECTRON  AFFINITIES 


Neutral 

H 

O 

F 

02 

03 

Li 


EA  (eV) 

0.75 

1.46 

3.40 

0.43 

2.10 

0.62 
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CALCULATIONS  ON  H30 


BASIS  SETS: 

4-31 G  (preliminary  survey) 
DZP  +  diffuse  s  and  p  on  O 


CASSCF  (Complete  Active  Space) 

9  active  electrons  (include  O  (2s) ) 

250  configurations  Cs  symmetry 

converge  on  1  2  A'  state 

solution  totally  symmetric 


FOCI  (First  Order  Cl) 

Single  excitations  with  respect  to  CAS 
DZP:  14902  configurations  2 A' 

SOCI  (Second  Order  Cl) 


DZP:  370330  configurations 


Calculations  to  date: 

molecule  in  C3v  symmetry 
calculations  in  Cs  symmetry 
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Vertical  Distance  d(a  ) 
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HgO  4-3 1G  MCSCF  Results  at  Fixed  H-H  Distance  ■  1.65a 
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INDICATIONS  FROM  PRELIMINARY 
4-31G  RESULTS 


1)  Can  locate  attractive  region  of  potential 
surface  due  to  ion-pair  bonding 

2)  Asymptotes  can  be  interpreted 

(  FOCI  -  a  =  1 .65  a0,  d  =  10.0  a0  ) 


-  h;  +  o-(2p) 

Hj  +  ot’D) 

2.25  eV  (spec  1.97) 
'  =—  H3  +  0(3P) 
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DZP/FOCI  RESULTS  AT  FIXED  H-H  DISTANCE 
AS  A  FUNCTION  OF  VERTICAL  DISTANCE 
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H3O  correlation  diagram 

PA'  Symmetry) 


H,0  LOWEST  POTENTIAL  SURFACE 
RELATIONSHIP  TO  OTHER  STUDIES 


(WD)  S.  P.  Watch  and  T.  H.  Dunning 

J.  Chem,  Phys.  TL  1303  (1980) 

POL-CI 

E.  Kraka  and  T.  H.  Dunning,  unpublished 

(NRS)  K.  S.  E.  Niblaeus,  B.  O.  Roos,  and  P.  E.  M.  Siegbahn 
Chem.  Phys.  25.  207  (1 977) 

UHF-SDCI 


barrier  to 


♦ 

6.1 

I 

oh*h2 


OHH2t 

(WD) 


Energies  (kcal m) 


dissociation 
to  H20  +  H 


* 


planar 

transition 

state 


(NRS) 
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present 
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GEOMETRY  DETERMINATION 
UHF  FOCI 


a  =  3.05ao  at  a  =  3.05aQ 
d  =  0.62ao  d  =  0.92ao 

H30  wave  function  qualitatively  same 
H3O*  ♦  diffuse  s  on  O 
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h3o  conclusions  to  date 

C3v  Restricted  Geometries 
DZP/FOCI 


Large  portions  of  potential  surface  may  be  characterized 
as  ion-pair 

Lowest  energy  (in  C3V  geometry)  corresponds  to  0(3s) 
Rydberg  and  dissociates  to  H20  +  H  with  small  barrier 


Ion-pair  minimum  (in  C3v  geometry )  at  ~a  =  3.05  a0, 
d  =  1 .0  aQ  106.5  kcal/m  above  H20  +  H 

No  adiabatic  correlation  with  ground  state 
OH  +  H2  or  H20  +  H 

Bound  with  respect  to  OH(A2I+)+  H2 
and  H20(A1B1)  +  H 

Second  ion-pair  local  minimum  at  a  =  1.65  a0,  d  =  3.5  a0 
143.9  kcal/m  above  H20  +  H 

Higher  in  energy  than  OH(A2l+)  +  H2 

Barrier  to  22A'  lower  minimum 


h3o  work  in  progress 


SOCI  for  better  determination  of  absolute  energies 


Characterization  of  surfaces 

At  most  C3v  geometries,  12A" 
equivalent  to  22A' 

12A"  correlates  to  ground  state  OH  +  H9 


frRTETS 

Characterization  of  surfaces 

Much  higher  in  energy  -  small  binding,  if 
any,  with  respect  to  H3  +  O  (3P) 


No  correlations  with  lowest  OH  +  H20  +  H 
asymptotes 
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High  Spin  States  of  CO  and  CH 
by 

Daniel  D.  Konowalow  and  Marcy  E.  Rosenkrantz 
State  University  of  New  York  at  Binghamton 
Binghamton,  N.Y.  13901 


Ab  Initio  Calculations  On  High  Spin  States  of  CO  and  CH 


by 

DANIEL  D.  KONOWALOW 
Department  of  Chemistry 
SUNY-Binghamton,  Binghamton,  NY  13901 

ABSTRACT 

As  part  of  o-ur  current  interest  In  potential  novel  energy 
a  to  rage  de  v  i  e«s-we  have  investigated  several  high  spin  states  of  CO 
and  CH.  The  results  of  our  calculations  on  the  -’L+  and  -’ll  states  of 
CO  and  the  **lT  state  of  CH  will  be  discussed.  The  possibilities  for 
spin-orbit  Interactions  involving  these  high-spin  states  will  also 


The  basic  idea  behind  this  phase  of  our  research  is  to  investigate  the 

properties  of  pairs  or  aggregates  of  high  spin  atoms  which  may  combine  to  form 

a  deeply  bound  low  spin  molecule.  The  binding  energy  of  the  ground  state 

molecule  may  in  principle  be  stored  as  suitably  protected  atoms  or  else  as 

3  + 

high-spin  van  der  Waals  molecules.  Examples  of  the  latter  are  the  2^  state 

of  H„  which  could  release  about  A. 5  eV  on  a  spin  flip  and  formation  of  the  + 

1  7  +  ° 

state  or  the  state  of  N ^  which  could  undergo  a  series  of  spin  flips  co 

form  the  ground  ^2  +  state  with  the  release  of  about  9.7  eV.  Figure  1  compares 

c  +  8  j  + 

our  1  E  van  der  Waals  curve  of  CO  together  with  the  ground  r  state 
potential  curve  to  emphasize  the  tremendous  energy  storage  potential  of  such  a 
high-spin  van  der  Waals  molecule. 

The  ^2*  state  of  ^  has  been  examined  from  this  point  of  view  by  Ferrante 
and  Stwalley  [J.  Chem.  Phys.  78,  3107  (1983)]  who  estimated  its  potential 
curve  by  adding  an  estimate  of  the  dispersion  attractive  energy  to  their 
calculated  repulsive  self  consistent  field  (SCF)  energy.  Their  SCF  plus 
dispersion  potential  had  a  of  just  over  40  cm  1 .  Later,  Partridge, 
Langhoff,  and  Bauschlicher  [J.  Chem.  Phys.  84,  6901  (1986)]  found  from  an 
exhaustive  set  of  configuration  interaction  (Cl)  computations  in  terms  of  a 
substantial  basis  set  and  careful  estimates  of  basis  set  superposition  errors 
that  the  D  value  was  in  fact  only  about  21  cm 

In  view  of  these  new  results  Ferrante  and  Stwalley  [J.  Chem.  Phys.  83 , 
1201  (1986)]  showed  that  their  earlier  values  of  critical  temperature  and 
triple  point  temperature,  for  example,  had  to  be  revised  downward  by  about  £ 
factor  of  two.  The  point  of  this  recitation  is  to  note  that  extreme  care  must 
be  taken  to  assure  the  highest  possible  accuracy  in  the  calculation  of  such 
weakly  bound  van  der  Waals  potential  curves  if  they  are  to  provide  useful 
estimates  of  thermal  properties,  for  example. 

We  are  currently  Investigating  high  spin  states  of  CO  which  promise  to 
give  greater  specific  impulse  than  high  spin  The  hulk  of  the  experimental 

and  theoretical  information  available  on  CO  treats  mainly  spin  singlets  and  to 
a  much  lesser  extent  spin  triplets.  The  only  published  information  on 
quintets  we  found  is  the  pioneering  work  of  O'Neil  and  Schaefer  (ONS)  fj. 


I 


(t_uid  4,01  )  ^§-t9ua'  Suiputg 


no 


Chom.  Phys.  S3,  3994  (1970)).  Figure  2  displays  .til  their  calculated  bound 


3  3 

states  corresponding  to  the  C(  P)  +  0(  P)  interaction  except  for  the  ground 


state.  Their  minimal  basis  set,  full  valence  configuration  interaction 


calculation  is  subject  to  very  substantial  basis  set  superposition  error  (we 


estimate  an  error  of  about  2S00  cm  '  at  R=4  a  for  the  1  state)  and,  thus, 

o 


can  provide  only  a  qualitative  guide  to  our  own  work.  We  show  in  Fig.  3  that 


the  ONS  curves  for  the  lowest  and  ^11  states  are  substantially  deeper  than 


our  own  (primarily  due  to  their  BSSE). 

Let  us  describe  our  own  calculations.  It  is  performed  with  the  apparatus 


shown  schematically  in  Fig.  4,  which  ably  is  operated  by  Marcy  Rosenkrantz  and 
Jim  Franci s-Bohr ,  among  others.  Wo  augmented  the  basis  set  of  Liu  and  McLean 
[unpublished]  with  two  s-,  two  p- ,  and  one  each  of  d-,f-  and  g-  functions 
needed  to  optimize  the  dipole,  quadrupole  and  oetupole  polar izabi  1  i te.s  of  the 
ground  state  atoms  C  and  0.  Our  polarizabi liLy  results  are  given  in  Fig.  5 
and  are  compared  with  available  literature  values  (in  parentheses)  [H.J. 

Werner  and  W.  Meyer  Phys .  Rev.  A1 3 ,  13  (1976);  K.A.  Reinsch  and  W.  Meyer  Phys. 
Rev.  A18,  1973  (1978)].  The  atomic  polarization  basis  is  needed  in  order  to 
insure  an  adequate  description  of  the  long-range  interaction  energies  commonly 


approximated  by  the  familiar  multipole  series: 


EM,p  ■  -C6R_6-CRK"B-C.0K'''"- 


Note  that  we  do  not  use  perturbation  theory  to  calculate  the  long-range 


interactions,  we  merely  use  its  familiar  language  in  this  discussion. 


We  have  found  that  there  is  a  substantial  amount  of  configuration  mixing 


between  the  two  lowest  }'  states  of  CO,  especialy  in  the  region  of  3-8  a  . 
Unless  we  were  extremely  careful  with  our  choice  of  reference  state(s)  and 
input  vectors  on  which  to  base  the  mu  1 t i ref erence  second  order  Cl  (SOCI) 


calculation,  we  obtained  potential  curves  that  appeared  to  be  nonsense.  Wo 
shall  not  regale  you  with  the  details  of  our  early  computations,  but  rather 


describe  the  approach  we  have  taken  to  obtain  our  current  best  results. 

r’  I  _ . .  *  .  . —  ii  — 1  r l i  '1  3  +  .  . —  i -  i. 


First  we  performed  SOCI  calculations  on  the  1,2  y.  states  where  the 
excitations  were  performed  from  the  seven  configurations  that  arise  from 


allowing  the  six  p  electrons  to  occupy  any  of  two  o  and  two  tt  orbitals 


consistent,  with  symmetry.  These  calculations  employed  the  molecular 


§ 


£ 


B 


SWS1 


W 


orbitals  obtained  from  an  SCF  computation  on  the  state  with  the 

3 

configuration  [  ]5a6alir  2ir  we  then  analyzed  the  natural  orbitals  as  shown  in 
Figs  6  and  7  which  suggested  that  a  good  multireference  SOCI  base  would  obtain 
from  a  CASMCSCF  in  which  the  six  "p"  electrons  were  allowed  to  occupy  any  of 
the  three  0  and  two  it  orbitals  consistent  with  symmetry.  The 
corresponding  SOCI  comprised  about  119,000  configuration  state  functions 
(CSF). 

The  resulting  binding  energies  for  the  1,2  and  "’ll  states  are  shown  in 
Fig.  8.  As  shown  in  Fig.  8,  we  obtain  slightly  different  results  depending  on 
whether  our  input  vectors  for  the  multireference  (119,000  CSF)  6/3  +2 
computation  came  from  the  SCF  computation  for  +  [  ]  5o  6o  1  tt  ^  2tt  or  from  a  state 
averaged  (1  and  2  6/3o+2?r  MCSCF  calculation.  We  presume  that  the  latter 

computation  is  the  more  reliable. 

Figure  9  shows  that  the  1  +  state  has  nine  or  more  bound,  J“0 

vibratonal  levels  and  that,  it  is  expected  to  be  relatively  stable  thermally  at 
or  below  room  temperature.  Kven  the  2  "’ll +  and  1  '’ll  states  are  seen  to  be  at 
least  weakly  bound  (the  former  has  at  least  two  bound  vibrational  levels),  so 
they  too  could  be  potential  energy  storage  states  at  sufficiently  low 
temperatures.  However,  it  remains  to  be  seen  what  the  radi at ive  lifetimes  of 
these  states  are. 

We  have  started  on  calculations  of  other  high  spin  quintet  and  triplet 

3  3 

states  ot  CO  which  correspond  to  the  C(  P)  +  0(  P)  asymptote.  When  tnose 
calculations  are  complete  we  shall  he  able  to  consider  various  decay 
mechanisms  of  the  potential  energy  s  orage  states. 


CM 


o 

o  o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

T— 

Y— 

CM 

ro 

IO 

1 

1 

1 

1 

I 

( I _ ujo)  /(6j8ua  6ujpujg 

.70 


Theoretical  Studies  of  Spin-Forbidden  Radiative 
Processes  and  Electronically  honadiabatic  Processes 
Using  ab  initio  Electronic  Structure  Methods 


James  0.  Jensen  and  Byron  H.  Lengsfield 
Ballistic  Research  Laboratory 
Aberdeen,  MD 


David  R.  Yarkony 

Department  of  Chemistry 
The  Johns  Hopkins  University 
Baltimore,  MD  21218 


Applications  of  a  recently  developed  methodology*  for  treating  spin- 


forbidden  radiative  processes  within  the  Breit-Pauli  approximation  will  be 


reported.  Briefly  the  electronic  wavef unctions  (Ij)  are  determined  through 
first  order  in  perturbation  theory,  Yj=  +  Y*  ,  with  S'®  and  Tj  expanded  in  a 


configuration  state  function  (CSF)  basis 


k  ck  ‘k 


1  \ 


C1  satisfies  the  secular  equation 


( H°-  E°)  C*  =0 


while  V 1  satisfies 


(H°-  L°)  V1 = 
^  1  ~ 


HSO  C* 


I 


The  CSF 


where  hso  is  the  full  mi  c  rosoop  i  e  spin-orbit  hami t t on i an . 
basis,  6,  is  developed  from  or t hunormal  molecular  orbitals  determined  within 
the  state  average  MCSCF  approximation. 

The  determination  of  using  lb  and  i  is  preferred  computationally  over 
the  traditional  eigenfunction  expansion  method  which  gives  as 


V 


1 


CKT  HsocI 

~  c1, 

,  T,  o  . .  o  ,  ~ 

iLi  -  Lk) 
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In  particular  use  of  eq.  3  permits  treatment  of  spin-forbidden  radiative 
processes  originating  in  coupling  to  bound  states  in  the  continuum 
('resonances').  The  use  of  eq .  4  in  this  instance  would  be  difficult  if  not 
imposs ibl u. 

The  situation  is  illustrated  with  the  determination  of  the  radiative 

1  3  —  —  2 

lifetime  for  the  a  A  >  XI  transition  in  CH  .  Okamura  et  al.  have  reported 

the  lifetime  for  this  transition  as  b.9  +  O.S,  -0.0s.  This  transition 

acquires  oscillator  strength  by  coupling  of  the  *  A  and  ^  T.  states  respectively 
3  1 

to  TT  and  F  states  embedded  in  the  CH  +  e  continuum.  The  solution  ol  eq.  _> 
will  be  obtained  in  a  large  CSF  basis  (~1U">  terms)  and  analyzed  using  a 
natural  orbital  procedure.  The  possibility  of  optimizing  molecular  orbitals 
to  describe  the  *  ’  *  ;T  'resonances'  using  an  iterative  natural  orbital  procedure 
will  be  discussed. 

Recently  the  radiative  lifetime  (?)  of  the  a^A  state  in  NCI  has  been 
measured  by  two  diiiermi  experimental  techniques  with  significantly  different 
results,  T=L44t)ms  and  ?--Zms.  To  address  this  discrepancy  calculations  of  the 
radiative  lifetimes  for  the  (  b  *  X  f  ,  nlA)  ►  X  Y.  transitions  in  NCI  have  been 


performed.  Our  results  support  the  longer  lifetime  measurement. 

If  time  permits  we  will  discuss  the  application  of  recently  developed 
computational  techniques  for  the  evaluation  of  first  derivative  nonadiabatic 
coupling  matrix  elements 


g(J,I,Ra,R) 


Vrf»r 

a  ~ 


using  analytic  gradient  methods  to  discuss  nonadiabatic  chemical  reactions. 
In  particular  we  will  consider  the  charge  exchange  reaction 


H+  + 


NO  -*■  H  +  NO 
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Spin  Forbidden  Radiative  Processes 

.  Change  total  spin  angular  momentum 
Breit-Pauli  Approximation 

.  Hso  spin  orbit  Hamiltonian 

I, so  =  hso  +  hsoo 

hso  le”  spin  orbit  interaction 
5oo  2e~  spin  other  arbit  interaction 
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See  "Spectra,  lifetimes,  and  kinetics  of  matrix-isolated  NI  b*E+  and  a^A  ", 
J  Chem  Phys,  84,  2907(1986). 


-NEXT  SLIDE- 


3  -  1 

See  "The  lifetimes  for  spontaneous  emission  from  the  X  L  (v=l)  and  a  A 
states  of  CH~" ,  J  Chem  Phys,  85,  1971(1986) 
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ADVANTAGES 

1.  Facilitates  incorporation  of  experimental  data 

2.  Nixed  rep  resent  at i ons 


Spaces  [SJH  A  DRY,  J.  Chem.  Phys.  83,  1168  (1985)] 


(H°  -  E j)  tj  *  -  Hs°*° 


Advantages 

1.  AVOIDS  (IMPOSSIBLE)  TASK  OF  OBTAINING  ALL 
ZEROTH  ORDER  EIGENSTATES 

2.  PERMITS  TREATMENT  OF  PROCESSES  ORIGINATING 

IN  COUPLING  TO  STATES  EMBEDDED  IN  A  CONTINUUM 


3.  ORBITAL  OPTIMIZATION! 


IN  GENERAL 


1.  TRANSITION  MOMENTS  ARE  SENSITIVE 
to  the  QUALITY  OP  THE 
WAVE  FUNCTION 


2.  SPACE  approach  permits  use  of 
a  large  CSF  space 


To  exploit  (2)  with  an  eye  to  (i)  we  have 

implemented 

SYMBOLIC  MATRIX  METHOD  for  H80  ip° 
DRY,  J.  Cheat.  Phys.  84,  2075  (1986) 


CONCEPT  : 

Representative  MATRIX  Elements 
DIRECT  PRODUCT  SPACES 

LIU  &  YOSHIMINE ,  J.  Cheat.  Phys.  74,  612  (1981 


Benefits 


AVOIDS  SIZE  LIMITATION  OF  CONVENTIONAL  Cl 
105  -  106  CSF's  'routine' 


PERMITS  MAXIMAL  EXPLOITATION  of  SPACE  METHOD 


;vj  VJT%j  wj  \rv  n  »V  FVWTIFKFTiTXrX. 
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CABLE  I.  Spin-lorbidden  transition  moments’  for  A  — .Vane  a  —  X  transition  in  NCI. 
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/X.,  id  1\X -2~)  fi,  (b  ‘l  +  ,Xyl~) 


•r CVB  [Ref.  4(a)  j  I  £*p 
U  Seeker  et  al.  (Ref.  5)  J 
^’avne-Colboum  (Ref.  6)  $Q<=- 
This  work  FOc  iOOO 
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Lengsfield  Jensen  Yarkony 
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GOAL: 


nff>  ir  la2  2cr2  3o  lw  3  VALENCE 

w  ♦  a  lo^2cr  3o^4o  lw 
na*3a  a*~ 4cr 

3  1  1 

Construct  *  n  space  appropriate  for  t 


STEP  1:  VARIATIONAL  PRINCIPLE 

3E“,  *A,  3*^n  into  State  averaged  MCSCA 

RESULT:  4o  is  very  diffuse 

CH(2n)  +  e 
as  expected 
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Step  2:  (a)  Construct  ,  a*A) 

♦  l ( 1 n,  x3^-) 


in  flexible  ‘’2nd  order  space” 

80,000  &  130,000  terms 

e 

(b)  Analyze  the  solutions  using 
NATURAL  ORBITALS 

QUAL I TATIVEL Y: 

?  ?  1 

1  a  2 °  3 ° 1 n  is  correct 


but 


la22a23a24al7T 


should  be  valence  orbital 


Step  3:  Using  4<*  NATURAL  ORBITAL  redo  Step  2. 


Conclusion:  Now  possible  to  optimize  the  space 

used  to  describe  the  perturbation 
in  a  manner  analogous  to  that  used 
to  optimize  space  for  the  zerotn 
order  wavefunction! 
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EXPERIMENTAL  STUDIES  ON  THE  SYNTHESIS  OF  NEW  NOBLE  GAS  FLUORIDES  AND  HIGH 
OXIDATION  STATE  ENERGETIC  FLUORINE  COMPOUNDS  INVOLVING  UNUSUAL  BONDING 
SITUATIONS 


W.  W.  WILSON  AND  K.  0.  CHRISTE 

ROCKETDYNE  DIVISION  OF  ROCKWELL  INTERNATIONAL,  CANOGA  PARK,  CA  91303 

The  primary  objective  of  this  program  is  to  demonstrate  the  feasibility 
of  synthesizing  hypervalent,  high  oxidation  state  fluorides  of  nitrogen, 
oxygen  and  noble  gases.  Typical  target  compounds  are  NF^,  OF^  and  ArF^. 

Our  initial  efforts  are  concentrated  on  the  synthesis  of  NF  ^ .  A  sapphire 
reactor  was  set  up  for  uv-photolysis  in  liquid  F^  at  -196°.  The  efficiency 
of  the  apparatus  was  demonstrated  by  preparing  NF^+BF^  from  NF^  +  F2  +  BF3* 
but  photolysis  of  NF^  +  F^  alone  did  not  produce  detectable  amounts  of  NF,_ 
or  NF^+F  .  A  matrix  isolation  apparatus  was  built  and  the  NF^-F^BF^ 
system  was  studied  in  Ar  matrix  at  5°K  using  uv-photolysis.  Since  no 
evidence  for  the  formation  of  either  any  new  species  or  NF^  BF^  was 
obtainable  under  these  conditions,  the  apparatus  was  modified  to  use 
microwave  discharge  instead  of  uv-photolysis  for  the  generation  of  F  atoms. 
This  method  has  the  advantages  of  producing  much  higher  F  atom  concentrations 
without  photolytically  destroying  the  desired  products.  The  efficiency  of 
the  system  was  demonstrated  for  the  02~F2  system  producing  the  known  02F 
radical  in  good  yields.  The  reactions  of  microwave  discharged  generated 
F  atoms  with  either  NF^  or  BF^  alone  in  Ar  matrix  so  far  have  not 
generated  detectable  amounts  of  any  new  species. 

18 

Preliminary  experiments  using  F  radio  tracer  studies  were  carried  out  to 

determine  whether  nitrogen  can  exhibit  a  coordination  number  of  five  with 

fluorine.  Unfortunately,  the  facilities  available  (at  the  University  of 

Glasgow)  did  not  permit  the  generation  of  sufficiently  large  amounts  of 
18 

the  F2  required  for  our  experiments. 


HYPERVALENCY  (Continued) 


NFc  COULD  EXIST  IN  COVALENT  OR  IONIC  FORM 
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NFr  CHEMISTRY.  EXPERIMENTAL  STUDIES 
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HEMISTRY.  EXPERIMENTAL  STUDIES 


POSSIBLE  ANSWER  FROM  10F  RADIO  TRACER  STUDIES 


NFr  CHEMISTRY.  F  RADIO  TRACER  STUDY 


PLANNED  WORK 


ROCKWELL  INTERNATIONAL 
ROCKETDVNE  DIVISION 
KM  Cwtfft  Aw./Cwoji  (Vk,  CA  91104 


1st  HIGH  ENERGY  DENSITY  MATTER  CONFERENCE 
Department  of  the  Air  Force  -  AF  Rocket  Propulsion 
May  12  and  13  1987,  Arlington,  Virginia- 

Contract 

NEW  HIGH  ENERGY  DENSITY  SMALL  RING  SYSTEMS 

Presentation 

THE  PREVALENCE  OF  RHOMBIC  STRUCTURES  IN  A2B2  TETRAATOMICS 

by  KOOP  LAMMERTSMA  AND  OSMAN  F.  GUNER 

Dept,  of  Chemistry,  University  of  Alabama  at  Birmingham, 
University  Station  219  PHS,  Birmingham,  AL  35294. 


l8t  HIGH  ENERGY  DENSITY  NATTER  CONFERENCE 

Department  of  the  Air  Force  -  AF  Rocket  Propulsion 
May  12  and  13  1987,  Arlington,  Virginia 

THE  PREVALENCE  OF  RHOMBIC  STRUCTURES  IN  *2B2  TBTRAATONICS 

by  KOOP  LAMMERTSMA  AND  OSMAN  F.  GUNER,  Dept,  of  Chemistry,  University  of 
Alabama  at  Birmingham,  University  Station  219  PRS,  Birmingham,  AL  35294. 

A  variety  of  reactive  solid  state  materials  contain  rhombic  building  blocks 

2- 

as  important  structural  ingredients.  Illustrative  are  the  formel  Al2Si2 
and  Li2x2"  (X  =  N,  P,  As)  fragments  in  the  Zintl'  complexes  CaAl2Si2  and 
MLi2X2  (M  =  Ce,  Zr,  Pr),  respectively,  whereas  the  formel  Be2Ge22-  part  in 
CaBe2Ge2  can  be  viewed  as  a  square  network.  These  building  blocks  are 
tetraatomic  structures  of  the  A2B2-type.  We  have  studied  by  ab  initio  MO 
theory  a  variety  of  10,  12,  14,  and  16  valence  electron  A2B2  species  to 
determine  whether  high  energy  density  rhombic  structures  are  prefered  at  the 
molecular  level. 

For  the  16  el.  species  C4,  C-jHB,  C2H2B2,  and  related  cations,  the  rhombic 
structures  are  indeed  minima  on  the  potential  energy  surface.  This  is  also 
the  case  for  the  14  el.  C3Be,  C3HBBe,  and  related  cations.  However,  for  the 
12  el.  C2Be2  the  linear  form  (triplet)  is  more  stable  than  the  rhombic 
isomer,  which  is  a  minimum  energy  structure. 

2  — 

In  our  full  potential  energy  search  on  the  limiting  12  el.  C2Li2  and  10 
el.  b2b  ;2  species  we  focus  on  the  prevalence  of  the  highly  energetic  rhombic 
structures  and  their  unique  bonding  properties.  Similarly,  preliminary  data 
on  the  16  el.  disilicondicarbide  C2si2  will  be  presented. 
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6-31G*  Geometry 
CJ-C.J  1.543  k 
C2-C2  1.502  A 


No  electron  density  between  and 

The  HOMO  (I)  is  a  non -bonding  <T-orbital 

The  short  distance  between  and  results 

from  two  CT-bridged  -Tf-bonding  orbitals  (II) 

I  H 
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C  - 


INVERTED  sp2  1+ 

H 


INVERTED  sp3 


cr  Deficient  bonding 


Important  Mo's 


Homo /a/ &  Or* 

-  2E  'rtHoriATtC 


V  nj».  r\*  ryjin*  rv".  t  w  v\:  v*  1  w-  - 

Table  1.  6-31G*  energies 

(in  tu),  heat a  of  formation  (ABfl 

1  kcal/aol),  and  atrain  energies  (SB,  in 

kcal/mol )  for 

four-membered  1 

ring  structures. 

compound 

formula 

energy 

ABf 

SE 

f  21 29 

C4 

-151.14598 

228 

1  2329 

C4H22+ 

-151.53998 

726 

Cyclobutadiene  dication  (26) 

35  C4H42+ 

-152.91650 

623 

Bicyclol 1.1.0) butene  (8)14 

C4H4 

-153.58016 

142 

127 

1,3-Cyclobutadiene  (15)^ 

C4H4 

-153.6412 

103 

65 

Bicyclo[1.1.0jbutane  (2)14 

C4H6 

-154.87176 

51.9 

65 

Cyclobutene  ( 27,)  ^ 

C4H6 

-154.89962 

37.5 

29  * 

Cyclobutane  (28)^9 

C4H8 

-156.09703 

6.8 

26 

14 

Hydrogen 

s 

h2 

-1.12683 

Enthalpies  given  to  the  nearest  tenth  are  experimental  data 
integer  values  are  based  on  calculated  energies.  The  strain  energy 
SE,  are  derived  from  Franklin  group  equivalents:  CB2‘£“4.93 
CB=0.80,  015-08=08=18.88,  and  C=C=24.57  kcal/mol. 
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Relative  (in  kcal/aol)  Energies  of  C<H22+,  C4H+, 
and  C4  Isomers . 
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compound  energy  geometry  overlap 

r(Ci-C3;  r(C-X) ,rCC-Y)  population 


1 1*™* 


CftVf 


4,MH,»Tr 


"it?*!? 


l-t  »  .  I.*  ***  A.*  im' 
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See  "Site  Preferences  and  Bond  Length  Differences  in  CaAl^S^-Type  Zintl 
Compounds",  J  Am  Chem  Soc,  108,  1876(1986). 
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See  reference  from  previous  slide,  but  page  1884. 
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PHOTOEXCITATION  OF  TRIATOMIC  HYDROGEN 


H.  Helm  and  L.  J.  Lembo 

Molecular  Physics  Department 
SRI  International,  Menlo  Park,  Ca  94025 

1.  Emission  Spectra  (Herzberg) 

2.  Formation  of  H3  in  Charge  Transfer 

3.  Photoionization  Spectra 

4.  Two  -  Photon  Ionization 

5.  Excited  -  State  Charge  Transfer 

6.  Future  Work 

Supported  by  the  Air  Force  Office  of  Scientific  Research 
under  Contract  FQ8671-8700  432 


ABSTRACT 


Submitted  to 

The  First  Annual  High  Energy  Density  Matter  Conference 
(12-13  May)  1987,  Rosslyn  Virginia 

PHOTOEXCITATION  OF  TRIATOMIC  HYDROGEN* 

H.  Helm,  and  L.  J.  Lembo 
Molecular  Physics  Department 
SRI  International  Menlo  Park,  Ca  94025. 

We  have  used  a  coaxial  laser  fast-neutral -beam  spectrometer  to 
study  photoabsorption  of  long-lived  species  of  Hj.  The  H-j  molecules  are 
formed  by  electron  transfer  from  cesium  to  mass-selected  H-j+.  The 
neutral  molecules  are  detected  by  photoionization  several  microseconds 
after  their  formation.  Single -photon  ionization  of  Hj  has  been  observed 
in  the  near  UV  (3620  to  3380  A)  via  autoionizing  Rydberg  states.  In 
addition  one -photon  resonant  two-photon  ionization  has  been  observed  via 
the  3d^A^'  state  at  5781  A.  The  photoionization  data  suggest  that  the 

O 

long  lived  H-j  species  are  molecules  in  the  previously  known  2p  A2" 
state.  We  will  discuss  these  results  and  our  current  understanding  of 
the  structure  of  this  triatomic  hydrogen  species. 

*  Work  supported  by  the  Air  Force  Office  of  Scientific  Research. 


Three  of  the  Emission  bands  observed  by  Herzberg  and  coworkers  are 
shown  by  the  full  arrows.  (Ref.  1*5) 

The  photoionization  step  from  the  metastable  H3  state  is  shown  by  the 
open  arrow.  (Ref.  6). 


Predissociation  of  n=2  states 


2s  vibrational  coupling 
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Predissociation  of  n=2  states 


Rotational  coupling 


Metastable  triatomlc  hydrogen 


Lifetime  of  2pA£  state  depends  on  the  degree  of  rotation  of  the  molecule 
(Herzberg,  Hougen  and  Watson,  Ref. 5) 
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The  residual  width  of  the  lowest  level  (N*0.  K=0) 
is  attributed  to  Doppler  Broadening 


Radiative  lifetime  for  2pA"  2-  2sA'  i  transition  (1090  cm  _1 ) 
has  been  calculated  to  be  b7  us  (Gellene  and  Porter,  Ref.7) 
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Formation  of  H  3  by  Charge  Transfer 


Ionization  Potential  of  Cs  :  3.89  eV 

Binding  Energy  of  n-2  Electron  to  H3+  ~  Ry/4 


References  7-13 


Photoionization  Experiment 
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Photoionization  spectrum  at  high  n*values 
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Autoionization  Mechanism: 

vibrationally  excited  core  and  bound  Rydberg  electron 
exchange  energy. 

The  peaks  B  and  D  (B’  and  D')  are  identified  as  transitions 
to  the  7d  (7s)  Rydberg  states  of  vibrationally  excited  H 
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Two-Photon  Ionization  of  H 
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ELECTRON  TRANSFER  FROM  EXCITED  TARGET 


ground  state  target 

Cs(6s)  +  H3  - ►  Cs++ H3  (n»2)  -  228  meV 

excitation 

Cs(6s)  +  hv  (3401  A)  - ►  Cs(1 1  p) 


example  of  transfer  from  excited  target 


Cs(11p)  +  H,+ 


Cs++  H3  (n=26)  -  228  meV 


detect  Rydberg  molecule  by  field  Ionization 


H3(n=26)+  £ 
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Future  work 


-  Identify  lowest  autoionizing  levels  in  the  Rydberg  series 

to  support  assignment  of  vibrational  states  of  the 
metastable  molecule 

-  Measure  the  linewidth  of  the  transitions  to  determine  the 

autoionization  lifetimes. 

-  Two-Color  Ionization  experiments  using  as  intermediate 

the  n=3  levels  observed  by  Herzberg  to  determine 
the  quantum  defects  of  the  np  and  nf  Rydberg  series. 

-  Search  for  metastable  species  of  the  Hs  molecule  for  which 

a  number  of  metastable  states  have  been  predicted 
theoretically  (Ref.  14). 
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